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Abstract 

We deal with the incompressible Navier-Stokes equations, in two and three dimensions, when 
some vortex patches are prescribed as initial data i.e. when there is an internal boundary across 
which the vorticity is discontinuous. We show -thanks to an asymptotic expansion- that there 
is a sharp but smooth variation of the fluid vorticity into a internal layer moving with the flow 
of the Euler equations; as long as this later exists and as t << where v is the viscosity 



1 Introduction and overview of the results 

In this paper we deal with the equations: 



d t v u + v" ■ Vv u + Vp" = vAv u (1) 
div v u = 0, (2) 

(N 

where the spatial derivative x runs into the whole space M d for d = 2 or 3, v v and p v respectively 
denote the velocity and the pressure of a fluid and v ^ is the viscosity coefficient. The derivative 
fields i>^-V and A are with respect to x and applied component- wise to the vector- valued function v u 
whereas 'Vp u denotes the gradient (also with respect to x) of the scalar function p v . The equations 
(OP)-© are respectively Newton's law of motion and mass conservation for a homogeneous (constant 
density) fluid so that ([2]) reads as an incompressibility condition. When v = the equations (HI)-© 
are referred as the Euler equations whereas v > corresponds to the Navier-Stokes equations. 

The study of these equations naturally depends on the function space in which initial data is 
provided, and in which solutions are sought. In this paper we deal with the class of initial data of 
the vortex patches. Roughly speaking such data correspond to Lipschitz initial velocity v\ whose 
vorticity luq := curl vo is discontinuous across a hypersurface of M. d , say given by the equation 
<Po(x) = 0. To be more precise we recall first the definition of Holder spaces: 

Definition 1.1 (Holder spaces). For an open subset O of M. d , for s in the set N of all natural 
numbers including and a real < r < 1 the Holder space C s ' r (0) of the functions of class C s {0) 
such that 

\d a u(x) — d a u(y)\ 

IMIc<m-(0) : = su P\ a \^ s {\\d a u\\ Lxl{0) + su Px ^ ye0 ; _ ) < +oo, 

\x y\ 

and Cc' r (0) the subset of the functions in C s ' r (0) which are compactly supported. 



*Laboratoire Jacques-Louis Lions Universite Pierre et Marie Curie - Paris 6; 175 Rue du Chevaleret 75013 Paris, 
FRANCE 



1 



We now define precisely what we mean in this paper by vortex patches as initial data. 

Definition 1.2 (Vortex patches). 1. Smoothness of the boundary. Let be given an open 
subset 0o.+ ofM. d of class C s+1 ' r . This means that there exists a function ipo G C s+1,r (R d , M) 
such that an equation of the boundary <90o,+ is given by <9C>o,+ = = 0}, such that 
Oo,± = {±¥>o > 0}, where Oq- is the interior of the complementary of Oq i+ , and Vipo ^ 
in a neighborhood of dOo t+ . 

2. Boundedness of the boundary. We assume that the boundary <90o,+ * s bounded. 

3. Smoothness of the velocity and of the vorticity. We consider an initial velocity vq in 
L^(R ) -the space of Lf oc (M. d ) divergence free vector fields- whose vorticity ujo := curl vq is 
in the Holder space Cc' r (0o,±), that is a compactly supported vorticity which is C s,r on each 
side of the internal boundary <90o,+ = {</?o = 0}. 

Such an initial configuration is referred as a vortex patch since in the two dimensional case 
d = 2, such an initial vorticity loq can be for example the characteristic function of a bounded 
domain of class C s+1,r , as it was introduced by Majda in [89J (however the more special case of the 
ellipses was already considered by Kirchoff). 

Remark 1.1. One should notice that in the 3d case the initial velocity vq is automatically in 
Ll(R d ) whereas in the 2d case it is in Ll(R d ) if and only if the total mass of the vorticity vanishes. 

Numerous results about existence of solutions for such data (and even for more general ones) 
are available either for the case of the Euler equations and for the Navier-Stokes equations. The 
main goal of this paper is to obtain an expansion for the solutions of the Navier-Stokes equations 
which describes as well as possible their behaviour with respect not only to the variables t, x but 
also to the viscosity coefficient v, in the limit vt — * 0. However it is useful to reexamine first the 
case of the Euler equations. 

1.1 The inviscid case 

We first gather from the literature the following compendium of results regarding the inviscid case: 

Theorem 1.1 (Chemin, Gamblin and X. Saint-Raymond, P. Zhang, C. Huang, P. Serfati ). For 
initial velocities vq as described in the definition \1.2\ the following holds true: 

1. Existence and uniqueness. There exists T > (which can be taken arbitrarily large when 
d = 2) and a unique solution v° G L°°(0, T; Lip(M. d )) to the Euler equations: 

d t v° + v° ■ Vv° = Vp°, (3) 
div v° = 0, (4) 

with vq as initial velocity. From now on we denote D the vectorfield: 

D := d t + v° • V. (5) 

2. Propagation of smoothness of the vorticity. Moreover for each t G (0, T) the vorticity 

to := curl x v°(t) (6) 

is Cc' r (0±(t)) where 0±(t) are respectively the transported by the flow o/Oo,± o,t time t that 
is 0±(t) := <-f°(i, Co, ±) where X° is the flow of particle trajectories defined by dt<¥°(t,x) = 
v°(t,X°(t,x)) with initial data X°(0,x) = x. 
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3. Propagation of smoothness of the boundary. For each t € (0, T) the boundary dO + (t) 
of the domain 0±(t) is C s+1,r . 

4- Propagation of smoothness of the level function. For each t E (0,T) the boundary 
dO+(t) is given by the equation dO+(t) = {(p°(t, .) = 0} ; where 

ip e L oo (0,r;C 1 ' r (M d ))nL oo (0,r ; c s+1 ' r '(O±(t))) 

verifies 

Dip = 0, (7) 
V°\t=o = <Po- (8) 

Moreover 0±(t) = {±ip°(t, .) > 0} and there exists r\ > such that for ^ t ^ T , and x such 
that \ip°(t,x)\ < 7] the vector n(t,x) := V x ip°(t,x) satisfies n(t,x) / 0. 

5. Incompressible Rankine-Hugoniot. For each t £ (0, T) the function (tv° .n)(t, .) is C 0,r 
on{\<p°(t,.)\ <v}- 

6. Smoothness in time. Finally the internal boundary d(D+(t) is analytic with respect to time 
and the restrictions on each side of the boundary of the flow X° are also analytic with respect 
to time with values in C s+1 ' r . 

Remark 1.2. We stress that the vector n is not an unit vector even if it is so at t = 0, since it is 
stretched when time proceeds according to the equation: 

Dn = -\Vv°).n. (9) 

Let us start a bunch of comments about Theorem 11.11 with the two-dimensional case. In this 
case the existence and uniqueness of a solution of the Euler equations in the more general case of 
an initial vorticity which is a bounded function with compact support was proved by V. I. Yudovich 
in 1963 cf. [83]. The corresponding velocity field is log-lipschitzian and admits a bicontinuous flow 
X . Since in the two-dimensional case the vorticity is simply transported by the flow: 

u;°(i,;r (t,x)) =io {x) 

this implies that in the case of a vortex patch as initial data the vorticity uj° at time t remains a 
vortex patch relative to a domain which is homeomorphic to the initial domain. However if one only 
use that the initial vorticity is a bounded function with compact support, one can only get that 
the smoothness of Yudovich's flow X is exponentially decreasing: it is in C exp ~ at where a depends 
on the initial vorticity (cf. for example [9], Theorem 7.26). Therefore in the case of a vortex patch, 
Yudovich's approach only provide that the boundary of a vortex patch is in c exp ~ at . The point is 
to know how the smoothness of the boundary of the patch really evolves. Numerical experiments 
of Zabusky [115] in 1979 suggested that singularities of the boundary of the patches would develop, 
presumably in finite time whereas the ones of Buttke [TB] suggested a loss of smoothness. In 1979 
Majda [89] theoretically studied the evolution of the boundary of piecewise constant vortex patches, 
by a contour dynamic approach, announcing local-in-time existence and conjecturing that there are 
smooth initial curves such that the curve becomes nonrectifiable in finite time. Constantin and Titi 
[35] studied a quadratic approximation of the equation governing the evolution of the boundary for 
which S. Alinhac [7J found some evidence of finite time breakdown. 

As a consequence the global-in-time persistence of the initial C s+l ' T smoothness of the boundary 
was very surprising when proved by Chemin [29] (see also the notes |26j . [24] . [25] and the survey 
[32] ). Actually Chemin's approach (which already yielded local-in-time existence in [27] . [23] ) is 
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based on a 2-local paradifferential analysis that is on the propagation of smoothness with respect to 
non-smooth vectorfields. The idea of using regularity properties with respect to a family of vector 
fields goes back to the work of Bony cf. [12] and to the work of Alinhac [3] and Chemin [22] in the 
case of non-smooth vector fields. The vortex patch problem falls in this latter case since it holds 
smoothness properties with respect to vector fields tangential to the boundary of the patch which 
moves with the velocity of the fluid. 

Remark 1.3. This approach even allows to deal with more general cases since it could apply for 
example to initial data which are irregular with respect not only to one hypersurface {ip°(t, .) = 0} 
but to the whole foliation of the hypersurfaces {tp°(t, .) = a} for a£l. 

In [11] Bertozzi and Constantin succeed to recover global-in-time persistence of the C s+1,r 
smoothness of the boundary in the special of constant vortex patches by the contour dynamic 
approach. 

The persistence of piecewise smoothness of the vorticity (Holder regularity up to the boundary) 
was proved later first by Depauw in the case s = [16], [H] and by Huang [79] in the general case 
s in N with a Lagrangian approach. 

Chemin's approach was extended to the three dimensional case d = 3 by Gamblin and X. Saint- 
Raymond in [58j , but it is only a short time result since in the 3d case the vorticity ui° is stretched 
along particle trajectories according to the formula 

u°(t,X°(t,x)) := D x X°(t,x).LU (x) = to (x) • V w X°(t,x) 

which "solves" the equation 

Du° = u° ■ Vv°, (10) 

UJ°\ t =0 = U)Q. (11) 

A rough estimate of the vorticity stretching is given by 

||w (t)||L» < e /<>°H«P \\u) G \\ L oo (12) 

Of course in some specific situations the vorticity stretching vanishes so that the existence and 
regularity results are in fact global in time i.e. T > can be taken arbitrarily large as in the 2d 
case. For instance a particular situation is the one of an axisymmetric initial velocity, as considered 
by Gamblin and Saint-Raymond (and others). We will not specifically consider this case here. 

Indeed the statement of the 3d case in Theorem II .11 is not strictly given in [57J; roughly speaking 
Gamblin and Saint-Raymond deal with tangential smoothness in the case s = 0. The tangential 
smoothness in the general case s £ N was alluded in the comment (ii) of section l.d of Gamblin 
and Saint-Raymond and rigorously proved by P. Zhang in the couple of papers |117j . [116] written 
in Chinese. The persistence of piecewise C 0,r smoothness (the case s = 0) was proved by Huang 
in [80j by mean of a Lagrangian approach (see also [50] section 3.1). We didn't found any proof 
of the persistence of higher order piecewise C s ' r smoothness in the literature. This is one reason 
which leads us to now give a few details of the whole proof of the 3d case, for sake of complete- 
ness, including a proof of the persistence of piecewise C s ' r smoothness in the Eulerian framework 
of Chemin, Gamblin and Saint-Raymond. Another reason is that we will use some extensions of 
the subsequent methods in the following sections regarding the case of the Navier-Stokes equations. 

Theorem 11.11 gathers the results we will need regarding the inviscid case but many other works 
have been done around the case strictly addressed in Theorem ll.il Let us mention some papers by 
Chemin and Danchin about the two-dimensional case for vortex patches with singular boundary. By 
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using the pseudo-locality of the Paley-Littlewood theory and transport equations by Log-Lipschitz 
velocities they show that if the initial boundary is regular (C°' r in Chemin's works [30], [31] chapter 
9, generalized into C s ' r in Danchin's paper [10]) apart from a closed subset, it remains regular for 
all time apart from the closed subset transported by the flow. When the singularity of the boundary 
is a cusp the corresponding velocity is Lipschitz what allows Danchin cf. [43 j . [42] . [39] to prove 
the global stability of the cusp with conservation of the order (actually he gives a global result of 
persistence of conormal regularity with respect to vector fields vanishing at a singular point, which 
generalizes the structure of a cusp). 

On the other hand several numerical simulations, based on an adaptive multi-scale algorithm 
using wavelet interpolation, are presented, showing that corners are unstable either immediately 
becoming a cusp in the case of an initial sharp corner or immediately becoming flat in the case of 
an initial obtuse corner. 

In [36] (see also [15]) N. Depauw addresses the case when the domain of the fluid is a bounded 
subset of R 2 with smooth boundary. When the boundary of the vortex patch is away from the 
domain boundary, it will remain so under the evolution, and Chemin's result remains valid without 
restriction. When the vortex patch is tangent to the material boundary (for transversal intersection 
there is a counterexample of Bahouri and Chemin) , the author proves that if the initial vortex patch 
is of class C 1,r there exists a unique solution in this class of vortex patches at least up to some time 
T* > 0. He also proves local-in-time existence and uniqueness for several mutually tangent vortex 
patches in R 2 . 

In [50] Dutrifoy proves local-in-time existence and uniqueness in the case when the domain of 
the fluid is a bounded subset of R 3 with smooth boundary, when the boundary of the vortex patch 
is away from the domain boundary, without restriction and when the vortex patch is tangent to the 
material boundary, under a technical condition. In this later case, his method allowed to complete 
the previously mentioned 2d local-in-time result by Depauw into a global one (with a slight loss of 
smoothness). 

Finally the 2d result by Depauw was recently extended to global in time results by Huang in 
[81] without loss of smoothness. 



1.2 The viscous case 

Let us now consider the Navier-Stokes equations. Numerous works (see for example [97], 
[59] . [20] ) show how the additional term uAv u regularizes the initial dataQ Still to the author's 
knowledge the smoothing effect has so far been analyzed in a global and isotropic way, and this 
is precisely one main goal of this paper to describe how the smoothing effect is (micro-)localized 
in the case where vortex patches are prescribed as initial data. In effect we will show a conormal 
smoothing of the initial vorticity discontinuity into a layer of width \pul around the hypersurface 
{(f°(t, .) = 0} where the discontinuity has been transported at the time t by the flow of the Euler 
equations. Hence the fluid vorticity 

uj u : = curl v v (13) 

depends -locally (innerly)- on an extra "fast" scale: - (cf. section I4TT1) and will be described 
by an expansion of the form 

^(t,x)^n(t,x,^^), (14) 



where the viscous profile £l(t,x,X) admits some limits when the "fast" variable X goes to the 
infinities: X — > ±oo. 



It is even crucial in existence theories both in the approach of Leray and in the one of Kato. 
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The idea to associate a viscous profile to an inviscid discontinuity seems to date back to Rankine 
|98j and is widely known when the discontinuity is a shock as for instance in compressible fluid 
mechanics (see the recent achievements by Gues, Metivier, Williams and Zumbrun in [69], |68j, 
[66], [67], [M])- However since they are characteristic and conservative the vortex patches are 
very different from the shocks of the compressible fluid mechanic (which are noncharacteristic and 
dissipative cf. for instance [96J). We therefore would like to precise that we borrow the words 
"viscous profile" to the setting of shocks profile but that our setting is quite different. For instance 
extra scales involved are not the same in the two cases. 

Still we hope that the approach developped here should be extended to some other setting 
where smoothing of characteristic conormal singularities occurs included the gradients jumps of 
compressible fluids, which are studied in the inviscid case in [6], [5], [3], [2], |102j . [95] . [9^] . in 
meteorology cf. |21j . the domain walls in ferromagnetism cf. [72J... 

These precisions done let us go back to our vortex patch problem and describe the construction 
of the viscous profile f2 involved in the expansion (I14p . We will look for a viscous profile £1 of the 
form 

n(t, x, X) = u°(t, x) + n(t, x, X), 
where 0(t, x,X) denotes a perturbation local with respect to the extra scale X so that 

lim n(t,x,X) = 0. (15) 

X— >±oo 

Hence the Navier-Stokes vorticities u) v (t, x) will be described by an expansion of the form 
ui u (t, x) ~ u>°(t, x) + O v (t, x) where u"(t, x) := U(t, x,^ 



vt 



The dependence of the perturbation oj u on v encodes the "conormal self-similarity" of the 

layer. Pragmatically the consequences of the condition (|15p on the profile Cl(t,x,X) at the level of 
the function ui u are threefold: 

1. For any (t, v) G (0, T) x R+, ui u (t, x) — > when <p° — ► ±oo. This was actually our motivation 
to impose the condition (|15p on the profile Q(i, x, X): it sounds natural that the viscous layer 
is confined to the neighbourhood of the hypersurface where the inviscid discontinuity occurs. 

2. For any t G (0, T), for any x G R rf \ dO+(t), u} u (t, x) —* when v — > + . This consequence is 
directly linked with another strong underlying motivation to this work that is the issue of the 
inviscid limit of the Navier-Stokes equation to the Euler ones. The "strength" of this inviscid 
limit (that is the functional space where it holds) does not depend only on the presence or not 
of material boundaries but also on the smoothness of the initial data. Basically the more the 
initial data is smooth the more the convergence is strong. For smooth data the Navier-Stokes 
solutions are regular perturbations of the corresponding Euler solutions in the inviscid limit 
(see for instance Swann [109] . Kato [85j, [86]) and converge say in any Holder spaces with a rate 
of order vt. We also quote here Masmoudi [91j for a slight improvement. At the other end in 
2D when vortex sheets are prescribed as initial data one only knows the weak L 2 convergence 
(cf. [II]). The vortex patches are an intermediary case, first studied in 2D by Constantin 
and Wu in [36]. In [l] Abidi and Danchin found the optimal rate in L°°(0, T; L 2 (M. d )), and 
recently [91] extends this result to the 3D case. We also refer to the papers of Hmidi [76J, 
|77] to the study of 2D vortex patches (including with singular boundary). In this direction 
the novelty here is again to describe locally what happens (actually the analysis performed 
here gives optimal estimates of convergence rates in any spaces as a simple byproduct). 
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3. For any (x,v) £ R d x E* , d^foa:) -> when i -> 0+. This yields that the Navier-Stokes 
vorticities uj v have the same initial value than the Euler one up . Let us mention here that 
the analysis can be simplified if on the contrary we allow ourselves to choose the initial data 
for the Navier-Stokes vorticities u u since there exists some well-prepared data for which the 
viscous smoothing is already taken into account (cf . section 14. 101) . 

We will argue (cf. section 14. 2p that the corresponding expansion of the velocity is of the form 

v u (t,x) ~ v°(t,x) + y/rtV(t,x, (16) 

Vi/t 



with V(t, x, X) satisfying 

V(t,x,X)->Q when X -> ±oo. (17) 

We will be led (cf. section |4~3"|) to consider for the profile V(t, x, X) the linear partial differential 
equation: 

LV = (18) 

where the differential operator L is given by 

L = S-t(D + A) (19) 

where D is the vectorfield in ([5]), and £ and A are some operators of respective order 2 and acting 
formally on functions V(t,x, X) as follows: 

£V := ad\V + ^-d x V - iv, (20) 

AV:=V.V x v°-2 iV - V * v0) - n n. (21) 

a 

Here a denotes a function in the space 

B := L°°([0, T], C°' r (R d )) n L°° (0, T; C s ' r (0± (*))) (22) 

such that 

inf a = c> (23) 

[0,T]xf d 

and such that a = \n\ 2 when \(p°\ < rj. 

We expect a continuous transition of the viscous fluid velocity v v and of the viscous fluid vorticity 
<jj v (these are the Rankine-Hugoniot conditions), instead of the discontinuity of the inviscid vorticity 
These continuity conditions would be translated into the following Dirichlet-Neumann type 
transmission conditions for the profile V(t, x, X) on the internal boundary {X = 0} (cf. section 
ET5J) : for any (t,x) G (0,T), xR d 

y = and = ^ ^, 24 

a 

where the brackets denote the jump across {X = 0} that is [f(t, x,X)] := f\x=o+ ~ f\x=o~ an d 
lu± are some well-chosen extensions of u \a ± (t)- 

The transmission conditions (|24p are normal for the operator E which is elliptic with respect 
to X thanks to the condition (|23l) . Actually we point out here that because of its unbounded 
coefficient X the operator £ does not strictly enter in the classical theory of elliptic operators (with 
t, x as parameter through the coefficient a). However we will see that it shares their main features, 
at least for our purposes. For instance omitting, to simplify, the dependence on t (which is here 
only parametric) we have the following result (cf. section 14. 4p : 
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Proposition 1.1. For any f G L 2 (R d ,fl'- 1 (R)) and g G L 2 (R d ) there is exactly one solution 
V G L 2 (M rf , of the equation SV = f with the transmission conditions [V] = and [dxV] = g 

across {X = 0}. 

One can verify that it is actually the kind of equation of the layer created in the very particular 
case of the stationary 2d circular vortex patches (see pQ). In this case because of the symmetry 
there are neither convection nor stretching, and the norm of the normal vector is conserved so that 
the profile equation is simply an ODE, whose solutions involve a Gaussian function. 

Of course the full equation (|18|) is much more intricate. Roughly speaking for t > the equation 
(I18p is hyperbolic in t, x and parabolic in t, X; but degenerates for t = precisely into the previous 
elliptic equation. However we will show that the equation (|18p . with the transmission conditions 
(j24"j) and the conditions (fTTj) at infinity are well-posed. We stress that since the hypersurface 
{t = 0} is characteristic for the operator L none initial condition at t = has to be prescribe for 
the equation (|18p . We will use here a L 2 setting, for two reasons: first in view to future extensions 
we want to give a claim hopefully robust. In particular it is well-known since [15] that in (multi- 
dimensional) compressible fluid mechanics the inviscid system should be tackled in L 2 -type spaces. 
This will to robustness is also the reason why we choose to put the emphasis on the velocity in 
this presentation, more than on the vorticity. The second reason for a L 2 setting is linked to the 
degeneracy at t = of the equation (fT8|) . which leads to the existence of parasite solutions. For 
instance if we look for solutions V not depending on X and neglecting the term involving A the 
equation (|18p simplifies into the Fuschian differential equation tdtV = — which admits an infinity 
of solutions i.e. V(t) = for C G M. However only one is in L 2 (0,T), corresponding to C = 0; 

and we expect that the scaling is enough relevant to have a solution with L 2 (0,T) smoothness, 
even in the case of the full equation (I18p . Let us give a precise statement: denoting E\ the space 

E 1 : = L 2 ((0,T) x R d ,H 1 (R)) 

we will prove: 

Theorem 1.2. For any f G E[, for any g G L 2 ((0,T) x M. d ) there exists exactly one solution 
V(t,x,X) G Ei of LV = f with the transmission conditions ([V], [dxV]) = (0,g) on T := (0, T) x 
M. d x {0}. In addition the function v^H^C^, -)H^ 2 (R d xR) * s continuous on (0,T). 

The equation LV = f is satisfied in the sense of distributions on both sides T>± := (0,T) x 
R d x R| of the hypersurface V := (0,T) x R d x {0}. Since V is in E 1 the jump [V] is in L 2 (r). 
The sense given to the jump of the derivative is actually a part of the problem. The idea is 

to give some sense by using the equation put in a weak form thanks to Green's formula. We will 
explain this in details in section [331 

In the case of the transmission conditions (|24l) the source terms are orthogonal to n. It is then 
possible to use the uniqueness part of the previous theorem to prove that the function V(t, x, X) ■ 
n(t, x) vanishes identically. This orthogonality condition is self-consistent with the incompressibility 
condition (see section I4.7P and with the linearity of the equation (see section I4.3.2P . 

We are now interested in the smoothness of the solution V given by Theorem 1 1.2 1 It is judicious 
to look again at the associated elliptic problem first. We will prove that the solution inherits the 
smoothness with respect to the usual variables t, x from the coefficients; and which are piecewise 
smooth with respect to the fast variable X. To be more precise, let us denote p — S(M) the space 
of the functions f(X) whose restrictions to the half- lines M± are in the Schwartz space of rapidly 
decreasing functions, and A the space (of the functions f(t,x,X)): 

A := L°° ((0, T), C°' r (R d ,p - «S(R))) n L°° (o, T; C s ' r {0±(t),p- 5(R)) V 

In section S3] we will prove: 
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Theorem 1.3. The solution V(t,x,X) of the equation SV = f with the transmission conditions 
is in A. 

The main idea of the proof is to use a spectral localization with respect to x, which is here 
parametric. The point is that this process is compatible with some classical elliptic arguments used 
to get smoothness with respect to X. 

We will be able to prove the same for the full equation (| 18j) : 



Theorem 1.4. The solution V(t,x,X) of the equation with the transmission conditions (%\) 
is in A. 

On the opposite to Theorem 11.21 we will use here the particular properties of the equation (|18|) 
throught the point [6] of Theorem 11.11 To explain this let us define for any Frechet space E of 
functions depending on t, x and possibly on X the space 

Ed := {/ eE/3C> 0/ (^y)fceN is bounded in E}. 

Thanks to the point [6] of Theorem 1 1.1 1 we will be able to construct (see section [4.3.5|) the extensions 
lu± of the vorticities and the function a in the space Ed- As a consequence we will actually prove 
in section fl~6l that V is even in Ad- 

The vorticity profile 0, in the expansion (|14p is then constructed as 

n(t,x,X) := u^(t,x) + n(t,x) Ad x V(t,x,X) for ±X > 0. (25) 

If piecewise smoothness of the initial data is sufficient it is possible to go on with the expansion 
with respect to vt of the solutions of the Navier-Stokes equations. At the extreme limit if the 
initial data is piecewise smooth on each side of the interface {(p° = 0} -that is if s = +oo- then it 
is possible to write a complete formal asymptotic expansion of the vorticity of the form: 

u v (t,x)=Y 'vV n j (t, x, * ^if) ) + 0(Vrt°°), (26) 

where the first profile 0° is the one previously constructed: 0° := $7. This construction will be 
achieved in section 14.81 

The stability of these expansions will be tackled in section [5j To describe the results we 
introduce the set T of the families (f v (t, x)) v of the smooth functions such that for any s' £ N, for 
any r' G (0, 1), the sequence 

is bounded. 

Theorem 1.5. There exists vq > such that for < v < vq for allk G N for any (t, x) 6 (0, T) xM rf 
w "(t, x) = Y vV W(t, x, y0( ^5 ) ) + Vrt k+1 u u R (27) 

a VVt 

with (uj v r )o <v< v in T. 

Let us point out that if it is well-known that for any v > there exists T u > and a unique 
solution v u E L°°(0, T u \ Lip(R d )) solution of the equations ©-© with vq as initial velocity, The- 
orem [L5] proves that the lifetime T u > can be extended such that T u > T. The lifetime of such 
expansions is the one of the solution of the Euler equation ("the ground state") which traps the 
main part of the nonlinearity of the problem. 
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2 Compendium on the Littlewood-Paley theory 

In this section we gather some usual results of the Littlewood-Paley theory that we will need. We 
begin with the case of the whole space for which we refer for example to the books [9] and 
for a much more detailed expository. 

2.1 Dyadic decomposition 

We first recall the existence of a smooth dyadic partition of unity: there exist two radial functions 
(p and x valued in the interval [0, 1], belonging respectively to 2^(5(0,4/3)) and to P(C(3/4, 8/3)) 
(where 5(0,4/3) and (7(3/4,8/3) denote respectively the ball 5(0,4/3) := {||£|| R 3 < 4/3} and the 
annulus C(3/4, 8/3) := {3/4 < ||£|| R 3 < 8/3}; and T>(U) denotes the space of the smooth functions 
whose support is a compact included in U) such that 

V£gR 3 , x(0 + E <^ 2 ~^) = 1 

\j - j'\ ^ 2 =>• supp (f(2~ j .) n supp f(2~ j '.) = 0, 
j ^ 1 =^=^ supp x(2" j .) n supp <p{2~ j .) = 0. 

The Fourier transform T is defined on the space of integrable functions / G L 1 (IR 3 ) by 

Tf := / e- x *f(x)dx 

and extended in an automorphism of the space iS'(M 3 ) of the tempered distributions, which is the 
dual of the Schwartz space 5(M 3 ) of rapidly decreasing functions. We will use the non-homogeneous 
Littlewood-Paley decomposition (in ^'(IR 3 )) 

id= £ A ;< 

where the so-called dyadic blocks Aj corresponds to the Fourier multipliers: 

A_! := x(D) and Aj := p(2~ j D) for j ^ 0, 

that is 

A_iu(x) := / h(y)u{x - y)dy and A jU (x) := 2 3j / h(2 j y)u(x - y)dy for j > 0, (28) 

where h := J 7 " 1 ^ and h := T~ x x- We also introduce the low frequency cut-off operator Sj := 
zJfe^j— l Ar 

2.2 Besov spaces 

We now recall the definition of the Besov spaces B^ q on the whole space M 3 which are, for A G R 
(the smoothness index), p,q G [l,+oo] (respectively the integral-exponent and the sum-exponent), 
some Banach spaces defined by 

B* q (R 3 ) := {/ G S'(R 3 )/ ||/b pVR3) := ||(2^||A J /|| iP(M 3 ) )^_ 1 || i , < oo}. 

These spaces do not depend of the choice of the dyadic partition above. When p = q = +oo the 
Besov spaces B^ q are simply the Holder- Zygmund spaces and in particular for A G M + — N they 
coincide with the C s ' r spaces of the introduction in the sense that B^ oo(M 3 ) = C s ' r (R 3 ) where s 
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is the entire part of A and r := A — s. It is worth mentioning the the space L°°(R ) is continuously 
embedded in the Besov space B® 



'00,00 



L°°(M 3 )^BSo )0O (R 3 ). (29) 
Conversely the spaces for A > the spaces -B^^M 3 ) are continuously embedded in the space 



L°°(R d ), and the L°° norm can also be estimated by the following logarithmic interpolation in- 
equality: 

< £(11/11*0 (R3) J/|| B * (R3) ), (30) 



where we define for a and b strictly positive L(a, b) := aln(e+ |) -which is notably increasing both 
with respect with a and b. We also recall the way Fourier multipliers act on Besov spaces: if / is 
a smooth function such that for any multi-index a there exists an integer m G N such that 

V£gM 3 , \d a f(0\ ^ C Q (l + |£ir- |a| , (31) 

then for all A G R and p,q G [l,+oo], the operator f(D) is continuous from Bp q to Bp~ m . In 

|2\ * 



particular introducing the function A(C) := (x(C) + ICI ) 3 > where x is m Cq°(R ) positive equal to 
1 near 0, and the corresponding Fourier multiplier A := X(D) we get an one-parameter group of 
elliptic operators A 4 , for t G M, continuous from Bp q (R 3 ) to -Bp+*(R 3 ). 

2.3 Bony's paraproduct 

When v and w are two Holder distributions, we denote by T v w Bony's paraproduct of w by v: 

T v w := Sj-ivAjW (32) 

for which we have the following tame estimates, for A 6 1, denoting [|.||^ the norm in the space 

B x 

00,00 

H^uHU ~ IMIi 00 IMU for A G R, (33) 
\\{v-T v )w\\x < \\v\\x\\w\\ L <*> for A > 0. (34) 

We will also use the following commutator estimate (cf. [9] Lemma 2.92): if / is a smooth func- 
tion homogeneous of degree m away from a neighborhood of then the commutator [T a , f(D)] := 
T a f(D) — f(D)T a between a paraproduct T a and the Fourier multiplier f(D) can be estimated for 
any AgR and for any r G (0, 1) by 

\\[T a , f(D)]u\\x~m+r ^ IkHrNU ( 35 ) 

We also mention the following useful estimate for commutators of the form Rj := [v° • V, Aj]f (cf. 
[9] Lemma 2.93): 

sup 2^' A \\RjWlp < \\Vv°\\ Bd/p nLoo ||/bA for0<A<l + d/p, (36) 

J ^ 1 

sup 2* x \\RA\lp < \\Vv°\\ L «,\\f\\ B x + ||V/||ip||Vu°||oA-i for < A. (37) 
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2.4 Transport estimates 

We recall some useful transport estimates (cf. [9] Theorem 3.11): for s = —1 or s = 

ll/Wllc-(M3) < (ll/|t=olb-(M3) + f \\{dtf + v° ■ Vf)(r)\\ c ^ m e- cv ^dr)e cv ^\ (38) 

J o 

where 

V(t) := [ \\Vv°\\ L oo m ds 
Jo 

and for any s ^ 1 

\\f(t)\\c^ m < (||/|f=o||c..r(RS)+ ril(5 t / + ^ -V/)(r)|| CSjr(M 3 )e -^-^«dr)e c ^W, (39) 

JO 

where 




We refer for example to [9] Theorem 3.11 for a proof of the estimates (I38p - (|39p by a spectral 
localization that is by applying the dyadic block Aj to the equation and then making use of an 
energy method, taking care of the commutators thanks to Bony's paraproduct. 

2.5 Besov spaces on Lipschitz domains 

When $7 is a Lipschitz domain one usually define the Besov spaces B pq (Q) (see for instance [110] ) as 
the set of restrictions of all elements of i?p 9 (IR 3 ) in the sense of 2?'($7), the space of the distributions 
on O. In other words the spaces Bp q (Q) consist of exactly those distributions / £ T>'(JY) which have 
extensions belonging to Bp q (K. 3 ). Endowed with the quotient space norms they become Banach 
spaces. 

We will make use of extension operators. We give the following nice general result. 

Theorem 2.1 (Rychkov |101j ). Let £1 be a Lipschitz domain in M n with a bounded boundary. There 
exists a so-called universal extension operator that is a linear operator ext such that 

ext : B s p>g (n) -+ B s P)g (R n ) 

is continuous for any s£M, 0<p^oo,0<q^oo and satisfies (ext u)\q = u. 

One can use this result to give intrinsic characterizations of Besov spaces on Lipschitz domains. 
Dispa [U] have shown in particular that -as in the case of the whole space- the Besov spaces 
B So,oo(ty sti11 coincide with the Holder spaces: for A <E K+ - N, B^ )0O (Q) = C s ' r (£l) where s is the 
entire part of A and r := A — s. 

We will also use the following result about characteristic functions as pointwise multipliers in 
Besov spaces. 

Theorem 2.2 (Frazier and Jawerth [51]). The characteristic function xn of a Lipshitz domain 
£1 C K n with a bounded boundary is a pointwise multiplier in B p q (that is the map u \— > xo, u is 
bounded from B p (W 1 ) into itself) iff 

max (( 1), n( 1)) < s < -. 

p p ' p 

Let us stress that this theorem says in particular that the characteristic function of a Lipshitz 
domain is a pointwise multiplier in the Holder spaces C ,r . 
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3 On the proof of Theorem 11.11 



In this section we give a unified proof of the 3d case of Theorem 11.11 in the Eulerian framework of 
Chemin, Gamblin and Saint-Raymond including the persistence of piecewise C s ' r smoothness. 



3.1 Existence from a-priori estimates 

We first recall that local existence and uniqueness of solutions of the Euler equations when the 
initial velocity field is assumed to be in C s+1,r (K 3 ) n L 2 (R 3 ) (that is without a hypersurface of 
discontinuity) was proved by Chemin in the paper |28j. Moreover Bahouri and Dehman in |10j 
prove that the criterion obtained by Beale, Kato and Majda (singularities will not form in the flow 
as long as there is no rapid accumulation of vorticity, i.e., as long as the integral, in time, of the 
L°°-norm of vorticity remains bounded) remains valid in this setting. Actually Gamblin and X. 
Saint-Raymond begin their paper [ST] by collecting these results so that we also to refer to [ST] 
Theorem 2.8 (uniqueness) and Theorem 2.9 (existence). We also refer to the book [9] Theorem 7.1 
and 7.20 for a complete expository about local existence and uniqueness of solutions of the Euler 
equations in supercritical Besov spaces. 

To construct solutions of Theorem 11.11 we proceed by regularization of the initial data taking 
limits thanks to some a-priori estimates -uniform with respect to the regularization parameter- on 
the smooth solutions given by the results above. We refer to [57] Proposition 5.2 for this passage 
to the limit, and focus now to the way to get the a-priori estimates. Actually we will precisely 
study in the following sections very carefully what happens in the passage to the limit when we 
regularize the Euler equations into slightly viscous Navier-Stokes equations. 



3.2 First smoothness property of the boundary 

As direct application of the estimate (|38|) with s = we have the following estimate for the function 
ip° which satisfy the transport equation ©-(JS]) 

ll¥>°(t)llc°."(R3) < lko|b"^(R3)e C7V( * ) . (40) 



3.3 Initial conormal vectorfields 

Let us introduce the vectorfields we will look for smoothness with respect to. We first observe that 
at the initial time t = the vectorfields 








' d 3 (p ' 




-d 2 <po 


-d 3 tp 


, wl : = 





, wg : = 




_ «9 2 ^ _ 




_-3i9?o_ 








are tangent to the foliation of hypersurfaces given by (fo. As in our case we deal with only one 
initial hypersurface {ifo = 0} of singularity we also add the vectorfields 





-<9i((l - x)a*). 



-<9 2 (xxi) 
«9i((l - x)x\) 




to the previous triplet, where \ is a smooth function compactly supported identically equal to 1 
when \if (t, .)| < rj. We therefore get a set of five divergence free vectorfields, in C s ' r (]R 3 ) and 



conormal to the boundary <90 +j o of the initial vortex patch: w l .n = for any 1 ^ i ^ 5 on <90_| 
Moreover this set is maximal in the sense that it satisfies 



[W 



inf sup 

R3 0<i,i<5 



\Wn X Wl\ > 0. 
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We now illustrate that the initial vorticity ojq has "good" derivatives in the directions that 
are conormal to the boundary of the patch. To do this we first set up the notion of conormal 
derivatives. For u £ L°°(R 3 ) we define the distributions 

u>o t u '■= div (u)q ® u) = dj(w l Q yU). 

i 

where we denote the components of the M 3 -valued vector field Wq by Wqj, for j = 1,2,3 and the 
spatial derivatives dj = d Xj . Since the w l are divergence free there holds w l f u = w l ■ Vu when 
this last quantity is defined. To express iterated conormal derivatives we will denote for k £ N 

(w }fu:=w^}(w^}(...(w$ k }u))...) 

where := (/3i, ...,/?&) € {1, ••,5} fc . We are going to prove 

Proposition 3.1. The k-th order conormal derivatives of the initial vorticity for any integer k 
between 1 and s, that is the vector fields (wot) f or P := {Pi, ■ ■■,Pk) S {1, •-,5}*, are in the space 
C-^ r {R 3 ). 

Proof. To do this we first use Rychkov's extension theorem 12. II to get from the restrictions coq\o 0± 
on each side Oo,± of the patch two functions u>o : ± both in C s,r (M 3 ). Since the initial vorticity 
can be decomposed thanks to these extensions and to the characteristic functions xo ± °f Co,± as 
co = u>o,+Xo . + +^0,- Xo - it is sufficient to show that the (itfof/ wo,±XOo ± are m C _1 ' r (^ 3 )- 
since the vector fields w l are tangent to the patch there holds (f^ot) w o,±XOo ± = Xo ± ( u, ot)^ c<J o,±- 
Moreover thanks to Theorem I2.2I the characteristic functions xoo + are some pointwise multi- 
pliers in the Holder spaces C~ 1,r so that the question of the C _1 ' r (lR 3 ) regularity of the (lOot)^o 
finally reduces to the one of the (toot) ^0 ± which can be simply got from the paraproduct rules 



3.4 Time-dependent conormal vectorfields 

Since the vorticity satisfy the transport-stretching equation ()10p we expect that the tangential 
smoothness is conserved. However the boundary of the patch moves with the flow so that the 
notion of tangential smoothness depends itself of the solution. Indeed we define time-dependent 
conormal vectorfields w l by the following formula which imitates the dynamic of the vorticity to 

w\t,x) := {w i -V x X°){t,{X°{t,.))- 1 {x)) (41) 

so that they verify the equations 

d t w l + v° ■ Vw i = w i ■ Vv°, (42) 

w%=o = wl (43) 

By a straightforward estimate we can show that for any t the set of the vector fields w l {t) 
remains maximal. Indeed (cf. [57] Corollary 4.3) they satisfy the estimate 

mt)]~ L i m < e ^/oll- ll^ ( M3) [Wo]-i (M3) . (44) 

Their divergence satisfy the equations 

<9 t div w i + v° • Vdiv w i = div w i ■ Vv° 

and since they are divergence free at t = the vector fields w l remain divergence free when time 
proceeds. 
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Let us denote for any circular permutation of (i,j,k) of (1,2,3) the vector u l := A w . It 
comes out that the vector n satisfies 

n(t,X°(t,x)) = £ u%X°(t,x)). (45) 

In particular the vectorfields w l (t) remain conormal to the moving boundary dO+(t) when time 
proceeds: w l .n = for any 1 ^ i ^ 5 on dO+(t). 

The conormal vectorfields satisfy the transport-stretching equation (|42h so that we should expect 
that their smoothness persists as long as we have good estimates of the velocity. In particular we 
need Lipschitz estimate as well as conormal and piecewise estimates of the velocity. Next section 
is devoted to the way to deduce such estimates from the ones on the vorticity. 



3.5 Static estimates 

Let us start with the Lipschitz estimate of the velocity. As the Marcinkiewicz-Calderon-Zygmund 
theorem fails to give a Lipschitz estimate of the velocity simply by ||u; ||^oo n ^2, we need to consider 
a larger norm in the r.h.s. But using normal derivatives (i.e. along n) would lead to estimates 
non-uniform with respect to the understood regularization parameter (the vorticity is expected to 
be only bounded not discontinuous across the patch boundary). The argument -dating back from 
Chemin for the 2d case (see [32] , [26] , [23] , [25] , [23] )- consists to use the ellipticity of the div — curl 
system to establish a Lipschitz estimate for the velocity field with respect to the L°° norm and 
to conormal derivatives only of the vorticity . Moreover one can manage to involve the conormal 
derivatives only through a log, that is widely known helpful in a subsequent Gronwall lemma since 
an inequality of the form y' ^ y logy does not lead to a blow-up. In the 3d case Gamblin and 
Saint-Raymond (cf. [57J Lemma 3.5) also succeed to use the ellipticity to get rid of the derivatives 
normal to the boundary proving: 

Proposition 3.2. There holds 

\\v°\\Lip(R3) ^ ll^°llL 2 nL°°(i; 3 ) + ^(II w0 |Il°°(r 3 ) ; ll^llc^''): (46) 

where we define for a and b strictly positive L is the function defined in £30)) and the conormal Ceo 
norm denotes: 

\\u\\ c 0,r 1= |H| L oo (K ,3) + [W]~i (R3) ^ (IKIl0,r|M| L oo (B >3) + \\w % j u\\-l >r ). 

Scheme of proof. Actually since can easily be estimated by 

II _,0 II <- II, ,0 1 1 

(thanks to an appropriated splitting into near and far field in the integral formula of the Biot-Savart 
law) the delicate problem is to estimate the ||<9if°||_L°o . Let us continue to reduce the problem by 
introducing a function x i n C^°(R 3 ) positive and equal to 1 near and x{D) the highly smoothing 
corresponding Fourier multiplier. Let A(£) := (x(C) + IC| 2 ) 3 an< 4 A := A(D) the corresponding 
elliptic Fourier multiplier. Since the low frequencies part can simply be estimated by 

H^-x^A-VlU,,. < [| V °[| £00 (47) 

we are left to estimate the high frequencies part 

11^(1 -x(tf)A-> IU«< II ^djdkLJ ^. (48) 
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By quadratic forms technics Gamblin and Saint-Raymond (cf. [57] Lemma 3.5 ) succeed to write 

djd k = Aa jik . + E didpb^u,;, (49) 

l,i,p 

where the functions aj jk and b l - l k are build by partition of unity from local expressions of the form 

(w m x w n )j(w m x w n ) k 

a 3,k = \—^rz — ' (so) 

\W m X W n \ z 

l ■ P l,i k (w m ,w n ) 

b j,k ~ | w m x w n\2 ' ( 51 ) 

where the P^ k (w m ,w n ) are homogeneous polynomials of w m , w n of degree 7 and w m x w n doesn't 
vanish, such that 

ll%Jlc^( R 3) < C(\\[W(t)] L i m . IKIM 
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We then rewrite the function A 2 djd k ui° in (j48]) as 



A^djdkto = h + f 2 where h := (1 - X (£)A~ 2 K fe u,° and / 2 := A" 2 ^ - Ao i( s>°. (52) 

Since ||/i||l°°(r3) ^ II 6< - ,0 ||l°°(r 3 ) we are left with the estimate of ||/2||l°°(r3). We then apply the 
logarithmic interpolation inequality: 

l|/2||i°°( R 3) < L(\\f2\\ B o coo (u3), ||/2||o,r)- (53) 

and observe that 

ll/allso (R8) < ||w || LO 



ll/2||o,r < (l|[W(t)]Zi( K a) E lhilo,,) 19 |k°ll c o. 



which allow to conclude to the estimate (|46H . 

□ 

In the same vein we have that conormal estimates on the vorticity imply conormal estimates 
on the velocity. More precisely there holds: 



IKt«°llo,r < H«' i tw ||-l,r+||« ||LipCR 8 )ll wi ||o,r. (54) 
Scheme of proof. To prove the estimate (|54"|) . we make use of Bony's paraproduct writing 

w i ti; = («; i -T t<I 0V« +T to4 Vu (55) 

and estimating the first term by ||w ||Li p (R3)||u; l ||o,r thanks to the paraproduct rule ([51}) . Then to es- 
timate the second term we split once again the velocity v° into low and high frequencies parts so that 
we get on one hand the term T w i f x(D)A~ 2 v° , which is easily estimated by H^HLiptRsjH^lkrj and 
on the other hand a combination, with coefficients valued in {—1, 0, 1}, of the terms T w iA~ 2 d k Vu . 
To control these last terms we commute the paraproduct T w % with the operator A~~ 2 d k V: 

T wl A- 2 d k Vu}° = [T w i,A~ 2 d k V]uj + A~ 2 d fc ti/ | oj° + A- 2 d k V(T^ - w^oj . (56) 

Now the first term in the r.h.s. of (|56|) can be estimated by 1 1 1 1 Li p (M 3 ) 1 1 ' lt;i 1 1 o,t- thanks to the 
commutator estimate (|35p with m = A = 0. The second one is estimated simply by \\w i f w°||_i ir 
and the third one by ||f ||Lip(R 3 ) ||w*||o,r thanks to the paraproduct rule (f34l) . □ 
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3.6 C°' r Dynamic estimates 

With the previous estimates on the velocity in hand we can now prove some persistence of smooth- 
ness of the vorticity. Actually Gamblin and Saint-Raymond (cf. [57] Proposition 4.1) prove the 
estimate 

\\w (t)\\ c o,r H Q < e ^ (ll^ollc-o.r H j-j—— — Tj J, (57) 



where V(i) denotes V(t) := J * ||v ||x,i P , the constant C depends only on r. All the norms above are 
related to the whole space M 3 . 

Scheme of proof. We apply the transport estimate (I38p with s = to the conormal vectorfields 
f = w % which satisfy the transport-stretching equation (|42p - (l43p and estimate 

\\dtw l + v° ■ Vu>*||cro,rnR8) = \\w l ■ VuHcTo.rms) (58) 



by ([HID- This Y ields 

\\w\t)\\ c o,r m < (Kllc°-( K 3) + / (||^t^ ||-l,r + ||^ || L ip(K3 ) || U ; i ||o 1 r)(T)e- G ' F MdT)e G ' F W. (59) 



' n| „ , _° II , _4-ll,.0|L. _.JL„<IL -V^-CVWj^CV® 



By using the equations (|10[) -(|lin of the vorticity lu° and the equations (|42p -(|43|) of the conormal 
vectorfields w l we next observe that the field w l f up involved in the inequality above satisfies the 
equation 

dtiw* t w°) + v° ■ V{w i f up) = w i t ■ Vu°), (60) 

Kfo; )| t= o=4t^o. (61) 

We apply the estimate ([3*5]) with s = — 1 to the field / = w % f w° 

IK f w°(t)||c-i.-(R3) < (114 t wo||c-i.r(R3) + / IK t (w° • V V )(r)|b-i„ (M 3 ) e- cy M ( ir)e cv 'W.(62) 

J o 

Let us recall that the initial data f cjoIIc-^CR 3 ) should be estimated by Proposition 13. 1[ Now 
to estimate the integral we first notice that 

w i f ■ Vv°) := div (w° 8> K f w )) + div (C O u°) (63) 

with C := div (u°®w ! — u; 4 ® w ). Then we make use of the paraproduct to commute the divergence 
operator in the first term of the r.h.s. of (|63p 



div (a; K f u; )) = [d k ,T wi ^ ]u° k + «9 fe ((w l j a; - r w i tw >g) , 
the sum over A; being understood. This gives -thanks to the estimate ([3~4"]) and ([35 



[|div(o; ®(ti; i tw ))[| £7 _x,r (R 3 ) < ||^t^°llo,r|k | 



Now to estimate the second term in the r.h.s. of (I63p we first notice that £ is divergence free so 
that 

||div(C® w°)||c-i.'-(R3) i$ H w °llLip(]R3)||Cllc-i.'- 
Since the tame estimates ([3"3"]) and ([3"""]) yield 



< [|«; i t« ||o,r + l|w°[|r<»|h"lk 
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we finally get -using one more time the estimate by (]54p - the following estimate of the conormal 
derivative of the stretching term 

|Kt(^°- Vv°)\\ c -i,r m < ||u || Lip(R 3 ) (||^U^ ||-l 1 r + ||^ ||i-(M3 ) |K||o,r)- (64) 
Now we observe that for any ^ s ^ t ^ T the vorticity lu° obeys the following L°° estimate: 

e V{s)-V{t) ^ M S )\\L~>m ^ e V(t)-V(s)_ (g5) 

II w (*)IIl°°(r 3 ) 



Combining the estimates (|62|) . (|64[) and (|65|) leads to 

IK t^°(*)llc-i.r( K 3) /Ikot^llc 



< ( " u ,, , v " ; (66) 



II w WIIl°°(k 3 ) ||wo||l°° 
+ f ll« (r)|| Lip(K 3 )( ||^(r)|| c c, (R3) + H^ t^Wllc-i.^) cy(r) cv(t) 



We add the estimates (f6"6"() and ([5U]) and conclude by a Gronwall argument. □ 

3.7 Gamblin and Saint-Raymond's result 

Putting together the estimates (|T2^ - (|4^ - (|4l)l) - (j57l) . Proposition 13. II and thanks to the relation ((431) 
we conclude that there exists T > and a unique solution v° G L°°(0, T; Lip(R d )) to the Euler 
equations such that for each t G (0,T) the boundary <90+(i) is C 1,r , and can be given by the 
equation dO+(t) = W°(t,.) = 0}, where cp° G L°°(0, T; C 1 ' r (IR d )) verifies The transport 

properties yield 0±(i) = {±yj°(t, .) > 0} and that there exists rj > such that for \<p°\ < rj, t ^ T, 
there holds Vi^° 7^ 0. We deduce the continuity of the normal component of the vorticity (u}°.n)(t, .) 
including across dO+{t) = .) = 0} from the diver gence free of the vorticity and its conormal 

regularity. 

3.8 Iterated tangential regularity 

In this section we briefly explain how to prove the propagation of the higher order smoothness- 
that is the C s+l,r smoothness- of the boundary dO+(t) (the point [3] of Theorem II .ip . This was 
done by P. Zhang in the couple of papers |117j . |116j written in Chinese, by adapting the method 
used by Chemin in the 2d case [27]. By local inversion there exists a parametrization of the initial 
boundary dO+$ given by a function 

Mn,T 2 )£C s+1 > r (S 1 xS\R 3 ). (67) 

Then for each t G (0, T) the boundary dO+(t) is given by the function 

x(t,Ti,r 2 ) := *(i,X (Ti,T 2 )) (68) 

and our goal in this section is to show that this function x is in L°° (0, T; (j s + 1 , r (g 1 x S 1 , R 3 )), what 
means that the boundary dO+(t) is L°°(0, T; C s+1 ' r ). 

Since for k = 1,2 the vector field d Tk xo is tangent to dO +t o at the point xo(ti,T2) and since 
the set of the vector fields {wq,Wq,Wq} is maximal on <90+,o there exists some functions 

a k ( Tl ,T 2 ) := (a k 1 (n,T 2 ),a k 2 (n,T2),a k 3 (T 1 ,T 2 )) G C^S 1 x S\R 3 ) 

such that 

3 

<9 Tfc Xo(Tl,T 2 ) =^2a k (T 1 ,T 2 )w t (^ (n,T 2 )). 
i=l 
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As a consequence by applying the derivatives := d^d^ to the function x in (|68l) for 
a := 02) G N x N with |q| := a± + 02 ^ s + 1 we get 

H / x 

1=1 /3e{l,2,3} ! 

with p := (ft, ...,#), («>ot)^ := (^t).-^!) and where the functions are in C '^ 1 x S\M). 
Going back to the definition of the w l in (|41[) we see that 

(toj f ^)(^xo(ri,r 2 )) = w i (i,x(t,ri,r 2 )) 
and by iteration that the iterated derivatives ^(u)ot)^<-f^ x o( T i) r 2)) can be transformed into the 



iterated derivatives y(w\ w@ l J (t, x(t, n, r 2 )), where we denote := fli-i). We are thus 

led to look for L°°(0, T; C°' r ) smoothness for the (urf)^ w\ 

Since each u/f commute with the material derivative dt + v • V the (w\)P w l satisfy 

d t (w]fw i + v° ■ V{w]fw' 1 = (u>t) /J VVu°, (69) 
((tat)^'^)| t =o = (wojfwl (70) 



Applying the transport estimate (|39|) with s = to the equations (j69j) - (|70|) yields for any 
t€(0,T) 

WW/w'mo, < ([Ifaot^jlkr + T ||M) /3 '^V V (T)|| 0ir e- cv/ Mdr)e cl/ W, (71) 



where V(i) := f \\Vv°\\i J oo(^3^ds is now under control. Moreover since their definition in section 
3.31 we know that the initial conormal vectorfields Wq are in C S,T . Since \(3'\ ^ s the paraproduct 
rules yield that the initial data (u>of)^ w l is actually in C 0,r . As a consequence the estimate (|7T1) 
now simply reads 



||HT u; 4 (t)||o, r < C(l + / ||Hf «,*V«°(r)||o, r dr). (72) 

JO 

Now we can estimate the r.h.s. of (|72p by following the approach of Gamblin and Saint-Raymond 
mentioned in section [3~5l proving on one hand that that iterated conormal regularity for the vorticity 
implies iterated conormal regularity for the velocity and on the other hand that iterated conormal 
regularity for the velocity is preserved when time proceeds. For the first step we refer to |116j p387 
which shows that for any time the C°' r norm of the (tyf) w lJ \7v° can be estimated by the C~ 1,r 
norm of the iterated conormal derivatives (w;f)^a; of the vorticity. More precisely when proceeding 
by iteration on \/3'\ s and understooding the quantities already known from previous steps we 
have 

[I WAlo.r £ IIMMl-l.r + \\(w}fw%, r . 

Now in order to prove the second step we notice that applying the derivatives (tof/ to the equation 
()10p and using again that the w l \ commute with the material derivative dt + v • V yields that the 
{w^Yu) satisfy 

d t {w]fuj° + v° ■ V{w\fuj = (w]fuj Vv , (73) 
((w]fuj )\ t= o = (wotfuo (74) 
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We use again the transport estimate (|39|) but this time with s = — 1, which yields for any t £ (0, T) 

\\(wifu°(t)\U, r ^ (||Kt)%)||-l,r + f \\{w]f^Vv\r)\U. r e- cv ^dT)e cv{t \ (75) 

Jo 

Let us recall that the estimate of the initial data is performed in Proposition 13. 11 The estimate of 
the r.h.s. of (j73[) is done in [116] p388. It can be seen as an extension in the setting of iterated 
conormal derivatives of the estimate (|64|) of the conormal derivatives of the stretching term and 
reads as follows 

||H)^W||_i !r < ||H)^°||-l,r + \\Wfw%, r . 
Plugging this in the inequality (|75p yields an inequality of the form 

IIMAAtJH-Lr < C i l + T \\H)^°\\-l,r + \\Hfw%, r ) 

JO 

so that Gronwall argument leads to 

||(^t)^°||-l,,< \\Wfw%,r 

Now the estimate (|7"2"j) reads 

||H/^(i)||o, r < C(l + f \\(wifw%, r (r)dT), 

Jo 

summations over (3' and i being understood. Hence applying a Gronwall argument to the previous 
estimate yields that the (iff)^ w % are in L°°(0,T; C°' r ). As a consequence for any a with \a\ ^ s + 1 
we get dfx is in L°°(0, T; C°' r ), which proves that the boundary dO+{t) is L°°(0, T; C s+1 ' r ). 

3.9 Piecewise transport estimates 

In order to prove the propagation of piecewise regularity (in the next section) we will use piecewise 
estimate for the transport equations of the form: 

dtf + v°-Vf = g, (76) 
f\t=o = fo, (77) 

where fo, v and g are assumed to be given, the two last ones being time dependent. To obtain 
these piecewise transport estimates we will need the following result about extensions of functions 
preserving the Lipschitz property: 

Theorem 3.1 (McShane [92], Whitney [112] ). Let 0, be a subset ofM n . There exists an extension 
operator that is a linear operator ext such that 

ext : Lip (0) -► Lip (R n ) 

is continuous and satisfies (ext u)\n = u for any u in Lip (Vt). 

Let us now claim our result of piecewise transport estimates: 

Proposition 3.3. There holds the following piecewise estimate: for k = — 1 or k = 

\\f(t)\\c^ {0±m < C(t)(\\f \\ ck ,r i0o±) + f\\g{T)\\ C Kr {0±{T)) e~ cv ^) dT ) e cv ± (t)^ (7g) 

J 
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where V±(t) := /„* ||f°||Lip (o±(s))ds and for any k ^ 1 

\\f(t)\\c^(o ± (t)) < C(t)(||/ || ofc .r (0o±) + ril&(r)||c*.r (0± ( T)) e- 0V ^.±Wdr)e cv *.-.±W, (79) 

JO 

where Vfc,r,±(i) := Jo ll t;0 llc* ft ' r (C' ± (s))^ s where the constant C(t) depend only on the sup for s running 
over (0,t) of the Lipschitz norm of the domain 0+(s). 

Proof. Actually we will prove an estimate on the side 0+ totally independent of the one on the 
side O- . To fix the idea let us consider the case of the estimate on + . We make use of Theorem 
12.11 to get from the restrictions /o|o 0+ and g\o + (t) °f the initial data /o and of the source term g 
some extensions _fo,+ an d g+ defined in the whole space M 3 , with -for any k — T 

ll/o,+ llc fc . r (R 3 ) ^ C|l/o|o ,+ llc fc ^(Oo,+)' ll9+llc fc '''(K 3 ) ^ C(t)\\g\o + (t))\\c k > r (G + (t))> ( 80 ) 

where the constant C and C(t) depend only on the Lipschitz norm of the domains Oo,+ and 0+(t). 

To extend the velocity field v° we distinguish two cases. When fe = -1 or fc = we make use 
of McShane- Whitney extension Theorem 13.11 to get a Lipschitz extension u9 defined in the whole 
space M 3 of the restriction v \o + (t)- When k ^ 1 we use again Rychkov extension Theorem 12.11 to 
get from the restriction v\a + (t) an extension defined in the whole space Mr, with 

Ib+llc^CR 3 ) < c ( t )\\ v °\o + (t))\\c k ' r (o+(t))^ ( 81 ) 
Then we apply the estimates (f38|) and ([39]) to the transport equation 

d t f + + v° + -Vf + = g + , (82) 

f+\t=o = fo,+- (83) 

This yields for k = — 1 or = : 

ll/+(*)llc^(R3) < (ll/o >+ ||c7^(R3) + r ||9 + (r)||^,-(R3 ) e- c ^W a !T)e^W, (84) 

JO 

where V+(i) := J ||f+||Lip (R3)^ s an d for any fc ^ 1 

\\f+(t)\\c^ m < {\\h+\\c^ m + f ||5 + (r)||c-(R3)e- c ^-+Wdr) e c ^-+W, (85) 

J o 

where Vfc, r ,+(t) := JJ ||«5-||c*.'-(RS)ds. 

Next we observe that since the domain 0+(t) is the transported by the flow of the initial domain 
C +i o, the transport equation ([76j) - ([77j) is well-posed in U te ( 0j r)W x without any boundary 

condition. By restriction of the transport equation (|82|) - (|83|) the function f+\o + (t) also satisfies the 
transport equation ([76|) -([77 |) in 0+,o- Now we can conclude that f+\o + (t) = f\o+{t) by uniqueness 
(let say in L°° since it is sufficient for our purpose. Of course it is possible to prove uniqueness in 
a much wider context, and without use of Fourier analysis, see the work of L. Ambrosio). Making 
use of the estimates (|80l) - (!8TT) and doing the same thing on the other side CL we get the desired 
estimates. □ 

3.10 Propagation of piecewise regularity 

In this section we investigate the propagation of piecewise regularity for the Euler equations. Indeed 
we finish to prove the points [21 SI and [5] of Theorem II .li For the commodity of the reader we recall 
here the decomposition ([52*]) : 

K- 2 d 3 d k oj° =h + f 2 where h := (1 - X (D)A- 2 )a j>k u° and f 2 := A" 2 ]T 8^/:^°. 

l,i,m 
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Let us denote by Q the collection of the w % and of the coefficients a,j : k and b l -\ which appear in the 

decomposition above. For any m E N* we introduce the set Op m (respectively the set Op m ) the 
collection of the operators T (resp. T) of the form 

Tf -.= n a il n 1 a h ...a in n m f (res P . 27 ■- (w ii ])n 1 (w i ^)n 2 ...(w i -])n m f) 

Our strategy is to prove by iteration for ^ k ^ s that 

For 1 < i < 5, w i £ L°°(0,r;C fc ' r (O ± (t))) and w i \ u° G L°°(0, T; C k ~ 1,r (0±(t))), (86) 

v° 6 L°°(0,T;C k+1 > r (O ± (t))), (87) 
For |a| < k, d a Vv° = £/ + ^ T «^° + A ~ 2 dj Yl ( 88 ) 

where S is infinitely smoothing, T m (respectively T m ) is in Op m (resp. Op m ). Before to begin the 
proof by iteration let us explain briefly the choice of such a formulation. Of course the estimate 
dSZD implies that w° is in L°°(0, T; C k > r (0±(t))). Still to prove the estimate §7$ we will use 
transport features of the vorticity and it will turn out that it is useful to get some extra static 
estimates to deduce piecewise estimates on Vt>° from piecewise estimates on uP . For an operator 
the property to act continuously between some spaces of piecewise regularities is usually referred 
as a transmission property cf. Boutet de Monvel and followers |13| . |14j . |62j . |100j . In particular 
the references above prove that the operator T : w ^ Vu satisfies the transmission property in 
the case of smooth boundaries. In the present case the smoothness of the boundary is limited 
and we need to estimate carefully how the smoothness of the boundary is involved in the constant 
of continuity of the operator T above in C 0,r (0±(t)). In particular as a byproduct what follows 
will extend the previous results giving that for an open subset + of class C s+1,r the operator 
T is bounded from C s ' r (0±) into itself. Actually the analysis would be complicated here by the 
time-dependent setting. The key point is the decomposition (|88|) which basically allows to write 
any derivatives of the velocity as the sum of a smooth term (the low frequencies part), a local 
term involving derivatives of the vorticity, what allows direct piecewise estimate and a term which 
despite nonlocal involves (iterated) conormal derivatives (up to commutators) which better behave 
(than normal derivatives) when estimated in the whole space. 

Let us now start the proof by iteration of (|86|) - ([87]) - ([88]) . Regarding the step k = we al- 
ready know that the assertions concerning the w % and the w l t a; -the estimates (|86p ~ are satisfied 
from the estimates (I57p (which actually even give that the w l and the w % f u>° are respectively in 
L°°(0,r;C ' r (R 3 )) and L°°(0, T; C _1,r (IR 3 ))) whereas the decomposition ([88]) reduces to the de- 
composition (|52|) proved in section 13.51 and recalled above. As a consequence to conclude the step 
k = it only remains to prove the estimate (|87p . As we have said above we are going to use 
transport features of the vorticity, applying the piecewise transport estimate (fTSj) in C 0,r (O±(t)) 
to the vorticity / = ur with the stretching term g = ll>°.Vv as right hand side (we recall that 
the vorticity satisfies the equation (|lUp ). Since the functions V± are under control by the previous 
sections we get an estimate of the following form: 

\\v°\\c°,rfo ± (t)) ^ 1 + / ||w .Vv ||cfo,r( 0± ( T ))dr. (89) 
Jo 

Now let us use the decomposition (f52|) to get some extra static estimates to deduce piecewise 
estimates on Vv° from piecewise estimates on a; . From the low/high frequencies splitting of section 
13.51 estimates (|171) - (ji8l) the question reduces to prove that the operators A -2 dj dk act continuously 
in C 0,r (O±)(t). Indeed -since the operator x(D)A~ 2 is smoothing- we get 

ll/lllco.^OiCt)) ~ ll a j,fc|lc°.''(O ± (t))-ll w0 |lc .''(O ± (i))- 
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Now we use the local representation of the coefficients aj : k given in (I50p to estimate the norms 
ll a j,fc|lc , °. T- by means of the norms ||w;*|| c o,r and [W] which have already been estimated (cf. section 
13.61 and estimate (|44|) ). In the same way than for the aj t k that is by using the local representation 
of the coefficients b A given in (I5ip we can estimate the norms \\b l - j,||c°> r by some quantities already 
under control. As a consequence we simply estimate ||/2||c o > r (0±(t)) by ||/2||c°' r (R 3 ) which reduces 
to the ones of \\w l f oj° ||c- 1 '' r (R 3 ) which was already controlled in section [3761 Estimate ([57]) . Finally 
we get 

and since we have the tame estimates 

II, ,0 V7„,0ii <■ II, ,0 1 1 IIV7»,0|l _i_ll,,0|| IIV7„,0|| 

\\uj .vv \\co,r(o±(t)) ^ ll w \\c°< r (o±(t))-\\ Vu ||l°° + ||w -II V-u ||c°' r (C' ± (t)) 
we infer the following piecewise estimate of the stretching term: 

II, ,0 V7„,0|| <■ ||, ,0 1 1 

\\UJ .Vv \\c°,r(o±(t)) ^ ll W \\C ' r (O±(t)) 

then by applying a Gronwall argument to (|89|) we conclude that the vorticity 

lo° 6 L°°(0,r;C7 ' r (O±(t))). 

Using once more the static estimates yields the estimate (]87p and the proof of the step k = is 
complete. 

At this point we can use that the normal component of the vorticity (uj°.n)(t, .) including across 
dO + (t) = {ip°(t, .) = 0} is continuous (cf. section |3~7|) to get the point E]of Theorem ll.il for each 
t G (0,T) the function (w°.n)(t, .) is C°' r on .)| < V }. 

Let us now assume that (|86l) - (|871) - (|88l) hold for ^ k ^ s — 1. We are going to prove that 
(f86|) - (f87]) - ([88j) hold at the order k + 1. To do so we first apply the piecewise transport estimate 
(|79p to the conormal vectorfields / = w l with as respective r.h.s. g = w l .X7v° (we recall that the 
conormal vectorfields w l satisfy the equations (pl2"|) - (j4*3|) ). Since the functions Vk+i )T ,± are under 
control by the previous step we get an estimate of the following form: 

\W\\ C k+i, r{0± n )) <l + I \\w\Vv \\ C k+i,rr o± r T)) dT. (90) 

J 

We consider the extensions (to the whole space) w± and v ^. given by Theorem 1 2 . 1 1 from the restric- 
tions w l \o ± and v°|o ± (t) on each side of the patch of w % and v°. We denote by the curl of v±. 
We have therefore 

||w*.Vw°|| C .fc + l i r( CI± ( t )) < ||w! t .Vfj : || C .fc + l,r( IR 3). 

Now that we are dealing with functions in the full space we can use the paraproduct. Indeed we 
can adapt the estimate (|54p into the following one (with only modifications of the indexes in the 
proof of flMl)): 

\\ w ± t«±Hfc+l,r ^ ||«4 t^±l|fc,r + ||«±||Lip(R3)lk±l|fe+l,r- 

Since we do not know if the universal extension operator given by Theorem 12.11 acts continuously 
in L°° we first use the embedding ([291) before to go back to piecewise estimates through Besov-type 
spaces: 

||«4 t «±||fe+l,r ^ lk± t<*>±IU,r + ll«±l|l,r||lW±|U+l,r 

< II™* t UJ °\\c k ' r (0±{t)) + ll u °llc 1 ^(0 ± (i))ll^lc* fc + 1 '''(e> ± (t)) ( 91 ) 
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and since IkllcWOiC*)) is already controled the estimate (j90l) becomes 

||w*|| C fc+l,r( 0± ( t )) < 1 + / \\w l \\ C k+l,r( 0± ( T )) + \\W 1 ] U \\ C k,r( O± ( T - ) )dT. (92) 



To estimate the conormal derivatives w l f ui° in C k,r (0±) we apply the transport estimate (|79|) 
to the equations ()60p - ()6ip . At initial time the conormal derivatives Wq f ujq is in C k,r (0±) since 
C s ' r (C-i-) is an algebra, and the amplification factor involves only Vk, r ,± so that we infer an estimate 
of the form: 

IK t W°|| Ck ,r (0±(t)) < 1 + T |k 4 t (W° t V«°)[| 0fcl r (0±(T)) . (93) 

•/II 

Now to estimate the integrals above we proceed as follows. We consider once again the extensions 
w l ± and v± so that 

IK t (w° t v« )[| okl r (o±) < IK t (4 1 v4)|| cfc , r(R3) . 

Now we recall the decomposition (I63p which for the extensions reads 

u4 t (<4 • V^) := div (^ ® («4 f ^±)) + div (C± ® 4) (94) 

with := div (8> lOj. ~~ w ± <8><4)- The first term of the r.h.s. of (JM|) can be estimated proceeding 
in the same way than in section 13.61 This yields 

|| div (w^ ® (to± tw±))|| fc.r(R8) < ||«4 t4llc"=+i.''(R 3 )H w ±ll i ° - 
Now for the second term we have the following tame estimates: 

1 1 div (C± ® llo fc ^(IR3) ~ 11^4 t ^± llc7fe^(R3) + ||u4||c*+l,r( K 3). 

Thus using the continuity properties of the extension operator the estimate (I93D now becomes 



|K fa; ||cfc,r( 0± ( t )) < 1+ / |K t w °llc fc ''-(o ± (r)) + \\ wl \\c k + 1 - r (o±M)- ( 95 ) 

JO 

Combining with the estimate (I92p and up to a Gronwall argument we conclude that the conor- 
mal vector fields w l and the conormal derivatives of the vorticity w' L f u>° are respectively in 
L°°(0,T;C k+1 ' r (O±(t))) and in L°°(0,T;C k ' r (O±(t))). Thus we have already proved that the 
estimate holds at the order + 1. 

We are now going to show that u is in L°°(0,T;C k+1,r (O±(t))). In order to prove this we 
apply the estimate (f79|) to get 

\\u \\ck+i,rr O± (t)) < 1+ / l|w Vu || C fc + i,r (o±{r) )dr. (96) 

•/ii 

Now we estimate the integral above thanks to piecewise tame estimates. Let us mention a technical 
point: despite the zero-th order estimate 

\\ab\\c°.r(o±(t)) ^ ll a llco.''(e> ± (t))l|fr|U 00 + IklkHHIco^o^)) (97) 

is a straightforward consequence of Definition 11.11 we do not know if the higher order estimates 
(k G N) 

\\ ab \\c k -T(0±(t)) ^ ll a llc fc '''(C'±(t))ll & IU 00 + ll a lli oo ll ft llc fc .''(0 ± (t)) ( 98 ) 
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are correct. One should think to use the extension Theorem 12.11 to go back to the full plane case 
but we do not know if the universal extension operator given by Theorem 12.11 acts continuously in 
L°° . Still we can deduce from Definition 11.11 and Leibniz rule the estimates 

\\ a b\\c k ' r (o±(t)) < \\ a \\c k ' r (o±(t))\\b\\c k ~ 1 ' r (o±(t)) + \\ a \\c k - l ' r {o±(t))\\b\\c k ' r (o±(t)) (") 

and these estimates would be sufficient here, since the lower derivatives would be always estimated 
in some previous steps of the iteration process. Here we get -still understanding any controled 
quantities- 

1 + ||Vf \\c>°+i,r(G±(t)) + ll^ \\c k + 1 > r (0±(i))' (100) 

It is now time to stare at Vt>°. Derivating the decomposition (|88p at order k and eliminating the 
double derivatives thanks to the identity (|49|) we get the decomposition (|88|) at order k + 1. Now let 
us remark that thanks to the local representation (I50p -(|5ip we can estimate the C k+1 > r (0±{t)) norm 
of the coefficient Q, by a power of the C k+1,r (0±(t)) norms of the w l which have been previously 
controlled. We also know from section 13.81 that the iterated conormal derivatives (w\ )^oj° are in 
L°°(0,T; C -1 ' r (0±(t))). As a consequence we get 

||Vv || C fc+l,r( C ) ± ( t )) < 1 + \\u \\ C k + l,r( O± ( t )y (101) 

Plugging this in the estimate (jlOOH and then the resulting estimate in the estimate (|96l) yields 

||w || C fc+l,r( C)± ( t )) < 1 + / ||w || C .fc + l,r( C | ± ( T ))(iT. 

J 

so that lo° is in L°°(0, T; C k+l,r {0±(t))). Using again the previous static estimate we infer that 
the estimate (]87p holds at order k + 1. The iteration can be done and the point [2] of Theorem ll.il 
is therefore proved. 

Using the relation (|45p we next infer that n is in L°°(0, T; C s ' r (0±(t))) which together with the 
tangential estimates give that tp is in L°°(0, T; C s+1,r (0±(t))) and finish to prove the point [U of 
Theorem 11.11 

3.11 Analyticity 

In this section we say a few words about the final statement regarding analyticity w.r.t. time (the 
point [6] of Theorem II. ip . Actually the study of smoothness of the boundary w.r.t. time is already 
done by Chemin in its pioneering work [27] . He incorporates the material field D into its conormal 
fields and carries out iterated conormal derivations. He follows his approach in [28] and prove 
the smoothness w.r.t. time for classical solutions (vq in C 1,r ) in any dimensions and for Yudovich 
solutions (with bounded vorticity) in two dimensions. This result was extended into analyticity 
w.r.t. time by P. Serfati [103J (see also its doctoral thesis) in the case of classical solutions and in the 
present case of vortex patches, thanks to Lagrangian methods (by considering the Euler equations 
as an ODE for the flow). Let us mention papers [SB], [55] which recover by a Eulerian approach 
the result by Serfati of analyticity w.r.t. time of classical solutions and prove Gevrey-3 smoothness 
in time of Yudovich solutions. Finally we mention that such results also hold in the case of solid 
boundaries: the paper [57] of T. Kato proves the C°° smoothness in time for classical solutions in a 
smooth bounded domain in any dimension, the paper [60] extends Kato's result to analyticity and 
also proves the analyticity of the motion of a body immersed into a perfect incompressible fluid 
with C 2,r initial velocity. 
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4 Internal transition layer profiles 



In this section we look for an expansion for the solutions of the Navier-Stokes equations which 
describes as well as possible their behaviour with respect not only to the variables t, x but also to 
the viscosity coefficient v. The sections 14. II to I3~3l give the heuristic of the derivation of the profile 
problem -the equations (|18|) -(j24"|)- mentionned in the introductory section 11.21 More precisely in 
section 14.11 we will identify the inner fast scale as ^ ■ This means that the initial discontinuity 

of the vorticity is smoothed out into a layer of size \fvt around the hypersurface {<p° = 0} where the 
inviscid discontinuity occurs. In section ET2"1 we pay attention to the expected order of amplitudes of 
the velocity and pressure profiles. We will see that it is natural to associate to a vorticity expansion 
of the form: 



Ul (t, X) ~ il(t, X, 



a velocity expansion of the form: 



v v {t,x) ~ v°(t,x) + VrtV(t,x/°( t,X 



and a pressure expansion of the form: 



p v (t,x) ~ p°(t,x) + utP(t,x, lf [L l 



vt 

In section l4T31 we look for the profile equations. We choose there to deal with the velocity formulation 
of the the Navier-Stokes equations. In section 14.41 we study the problem obtained when setting 
formally t = in the profile problem. In section 14.51 we prove the existence and uniqueness of the 
layer profile in a L 2 setting, and we prove the smoothness properties of this solution in section 14.61 
The other properties of the profile are proved in section 14.71 what achieves the proof of Theorem 
11.41 given in introduction, n section 14.81 we write a complete asymptotic expansion. In section 
14.101 (respectively in section [4.11|) we mention how the analysis is simplified in a setting of "well- 
prepared" data (resp. in the setting of conormal singularities weaker than vortex patches). In 
section 14.121 we link the present study with previous works on stronger singularities. 



4.1 The inner fast scale 

The goal of this section is to explain how we identify the inner fast scale as 



■Jut ' 



4.1.1 A highly simplified model 

To identify the inner fast scale we first look at the Id scalar heat equation: dtuj u = vd 2 £ u v which 
plays here the role of "baby model" for the NS equations. We prescribe as initial data a discontin- 
uous vorticity: uf\t=Q = 1m + - In the inviscid case v = -which stands for (highly) simplified Euler 
equations- the solution is simply equal to the initial data ui°(t, .) := 1r + for any time, whereas for 
v > and t > one can explicitly compute the solutions oo u by convolution: 

r 1 r°° 2 

u u (t,x) := n(-=) where Q(X) := —= \ e~ y dy. (102) 



Hence the initial discontinuity of the vorticity is smoothed out into a layer of size v vt where occurs 
-smoothly- the transition between the values and 1. It is useful to rewrite the uj u as: 

(jj"(t,x) := u°(t,x) +&±(—=) when ± x > 0, 

\Jvt 
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where 



n±(X):=— e~ y2 dy when ± X > 0. 



One then see the "viscous" solutions u> u as the sum of the "inviscid" solution uj° plus a "double 
initial- (internal) boundary layer" Q± which satisfies the double-ODE: 

dx&± + y d x &± = when ± X > 0, (103) 
match the continuity conditions of oj u and c^uA at the internal boundary x = X = (for f > 0): 
^1^0+ + n + \ x=0+ = 1-1/2 = + 1/2 = u°\ x=0 - + fi_| x=0 - , (104) 

d x &+\x=o+ = 7Ti= = dx&-\x=o- (105) 

and vanish 

O±(X)^0 in the limits X ±oo. (106) 

These last limits correspond both to the limits i > 0, x — > ±oo and the limits ±x > 0, vt — > + 
(we recall that X is the placeholder for -^=)- Conversely the two second order elliptic equations 
(|103f) with the four "normal" boundary conditions (|104l )- (ll05p - (1106p (this last one contain two 
conditions) determine uniquely the profiles Cl± . 

4.1.2 Inner scale in the NS equations 

Of course the case of the NS equations is really much more complicated than the previous baby 
model. In particular the inviscid discontinuity moves: Theorem 11.11 states that the inviscid discon- 
tinuity occurs at the hypersurface {tp (t, .) = 0} given by the eikonal equations ©-(JH]) associated 
to the the particle derivative dt + v° ■ V x . Therefore we are led to consider the inner fast scale 
- and we expect that in the case of vortex patches as initial data the solutions oo u of NS can 
be described by an expansion of the form 

LO v (t,x) ~ L0°(t,x) + n(t,x, (107) 

yut 

where Cl(t, x,X) denotes a local perturbation so that 

lim Cl(t,x,X) = 0. (108) 

X^±oo 

Actually we will consider profiles even rapidly decreasing at infinity. 

Remark 4.1. At this point experts in mathematical geometric optics should argue that the fast 
scale y should be more intuitive since applying the particle derivative dt + v • to a function 

£l(t,x, - ) produces the singular term 



1 



which fortunately vanishes thanks to the eikonal equation. 
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The point is that such a choice of inner scale does not lead to totally satisfactory results. Let 
us show that on our baby model. If one look for representation of the solutions tu u of the form 

u"(t,x) := oj°(t,x) +n±(t, 4=) when ±£>0, (109) 

one then see that the "double (internal) boundary layer" fl± has this time to satisfy the parabolic 
equation: 

d t fL± = d 2 x Q. ± when ± X > 0, (110) 
with the same boundary conditions: 

n + \ x=0+ - n_\ x=0 - = -l, (in) 
d x n+\ x=0+ - d x n„\ x=0 - = o, (ii2) 

0±(-X*) -► in the limits X -> ±oo, (113) 

Since we prescribe the same initial data for oj u than for u , one has to prescribe zero initial data 
for the layers: 

n±|t=o = o, (ii4) 

so that the compatibility condition between the transmission condition (jllip and the initial con- 
dition (|114p on the "corner" {t = X = 0} is not satisfied even at order zero, which ruins any hope 
for smoothness regarding the profiles 

Now let us stress that if the "phase" ip(t,x) := v chosen in the asymptotic expansions 
(|107p does not satisfy exactly the eikonal equation it however satisfies 

(d t + v° -V x )^ = -^ (115) 

so that applying the particle derivative dt + v ■ V x to a function 

w w (t,x):=n(t,a;,^%^) (116) 



produces the term 



■±=(dt + v .v x )fi>.(d x n)\ v o M =-hxd x n)\ , 0{M) . (117) 

\ V x - ~m It A - 7775 



which is not singular (with respect to v) anymore. The point is that in the present setting of a 
localized profile Q the derivative Xd x in (|117|) does not cause any difficulty (note that such a term 
even appears in our baby model see (jl03jl ). in particular because the prefactor \ echoes the one in 
the term 

lniL ' r)l2 -d x n\^, x) (us) 



t a x = . 

which is the larger one produced by applying the Laplace operator vA x to the function u"(t, x) in 

(hud. 

Of course this additional derivative Xd x is more problematic in the traditional context of 
periodic oscillations of geometric optics. 
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4.2 Amplitudes 

We now pay attention to the expected order of amplitudes of velocity and pressure profiles. In the 
full plane the Biot-Savart law has Fourier symbol — A . It is a pseudo- local operator of order —1 
so that we expect that the velocity v u given by the Navier-Stokes equations can be described by 
an asymptotic expansion of the form: 



v"(t,x) ~ v"{t,x) + VvtV(t,x, v ' - ), (119) 

\/vt 



where the profile V(t,x,X) is also expected to satisfy 

lim V(t,x,X) = 0. (120) 

X— >±oo 

Arguably since the Euler velocity v is Lipschitz we expect that its viscous perturbations v u are 
uniformly Lipschitz, what gives support to the ansatz (|119p . Plugging (|107p and (|119p into the 
relations (I13p . taking into account ([6j) and equalling the leading order terms leads to 



nAd x V = n. (121) 

Hence the vorticity profile has to satisfy the orthogonality condition: 

n.n = 0. (122) 

This condition is not a surprise: since w° is divergence free w°.n is continuous (cf. the point [5] 
of Theorem ll.ip ) so that no (large amplitude) layer is expected on the normal component of the 
vorticity. 

Now the pressure p v can be recovered from the velocity v v by applying the operator divergence 
to the equation (pQ) which yields the Laplace problem: 

A, f = -div'.drf. (123) 

If the velocity v v satisfies the expansion (j!19|) . the r.h.s. of f)123|) should admit an expansion of the 
form: 

, &n _ f ^^) 



A xP u ~ A xP u + F(t,x, 

\Jvt 

where the function F satisfies 

lim F(t,x,X) = 0. 

X^±oo 

Since the Laplacian is of order —2 we are lead to consider a perturbation of order vt on the pressure: 

p v {t,x) ~ p°(t,x)+utP(t,x,^^-), (124) 

where -once again- the fast scale ^ is expected to be a local inner scale (the Laplacian operator 
is pseudo-local) so that 

lim P(t,x,X) = 0. (125) 

X^±oo 
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4.3 Looking for a profile problem 



Now that we have an intuition of the amplitude of the profiles, we look for the profile equations. 
We choose here to deal with the velocity formulation of the NS equations, which is believed a 
more robust method (in view of future adaptation to the compressible case for instance). We 
proceed in several steps. In section 14.3.11 we plug the ansatz into the velocity equation equalling 
the leading order terms. We then pay attention to the divergence free condition what leads to a 
crucial observation in section 14.3.21 In section 14.3.31 we get rid of the pressure in the velocity profile 
equation. As the vectorfield n may vanish, away the patch boundary, it is useful to modify the 
resulting equation in order to avoid a degeneracy of the order. This will be done in section 14.3.41 
In section 14.3.51 we study the transmission conditions on the inner interface X = 0. We will use 
several times the formulas (for any h) 



t 



(z4)2 h(t,x, 



<yt)l (dMt^^) + ^d x h(t, X ,^4) (126) 



ut 



ut 



ut 



V, 



h[t,x, r =-) 



v ,., <AVe), n(t,x) ( <P°(t,x) ^ 
V x h{t,x, 7 =—) H j^dxhfax, — 



ut 



(127) 



and 



uA n 



h(t, x, 



<p°(t,x) 



\n(t,x)\' 



d x h(t, x, 



<fP(t,x) 



+ ^_( A / d x h[t,x, 

+^ A x h(t,x,^l). 
t Wut 



ut ' 
<fP(t,x) 



(128) 



) + 2n(t, x) ■ X7 x dxh(t, x, 



<P°(t,x) 



)) 



4.3.1 Velocity equation 

We plug the ansatz (|119|) and p24|) into the equation (fT]), equalling the leading order terms, which 
are of order y/vf : 

dtv° + v° ■ Vv° + Vp° + (d t + v° ■ V x )v°.d x V = 0. (129) 

which is satisfied since the velocity v° satisfies the Euler equations and tp° satisfies the eikonal 
equation . At the following order \Jvi we get the equality 



3 t V + v° ■ V X V + V ■ nd x V + V ■ V x v° + d x Pn = ^{\n\ 2 d 2 x V + ^d x V - ^V). (130) 



4.3.2 Incompressible transparency 

We now pay attention to the divergence free condition. Plugging the ansatz (I119p into the equation 
([2]), retaining the terms at order y/ui and taking into account that the velocity v° given by Euler 
is divergence free leads to the orthogonality equation: n.d x V = 0, which by integration, with the 
condition (I120p leads to the condition: 

n.V = 0. (131) 

An important consequence of the condition (|13ip is to kill the third term in (|130f) which is the only 
nonlinear one. Here lies an analogy with the WKB theory of the propagation of high frequency 
oscillations for hyperbolic systems (see for instance |96j). The condition (I131D can be seen as a 
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polarization of the singularity on the components tangential to the "phase" ip . Then the vanishing 
of the burgers-like term in (|130p can be interpreted as a transparency property: the self-interaction 
of the singularity vanishes because this latter is characteristic for a field which is linearly degenerate 
(actually this concept belongs to the hyperbolic theory but the incompressible limit is reminiscent 
of this fact). 



4.3.3 Getting rid of the pressure 

The equation (|13U|) involves both V and P. However the pressure in the NS equations is not truly 
an unknown but can be recovered from the velocity (as recalled in (|123p ) so that we expect that the 
same holds for the profiles. One way to proceed is to project normally the equation (|13U|) . taking 
into account the equation ([9]) for n(t, x) and using the condition (I131|) . to get 

a x P:=-^- V f>- n . (132) 
\n\ z 

For ±X > 0, we integrate between X and ±oo, taking the condition at infinity (|125p into 
account, to find 

(V ■ V x v°).n , x 

P=-2 ± x ' . (133) 



jx l"l 

Remark 4.2. One another way to proceed is to explicit the term F which occurred in the expansion 
(|124p when we were discussing the amplitude of the pressure layers by mean of the Laplace problem 
satisfied by the viscous pressure p v . Taking into account the condition (|13ip . we get 

F = -2{d x V ■ V x v°).n. (134) 

Plugging on the other hand the expansion (|124p in the l.h.s. of (|124p yields 

\n\ 2 d x P = -2(d x V • V x v°).n (135) 

and we recover (|133p by integrating twice with vanishing conditions at infinity for P and d x P. 

It is interesting to note that despite the second method involves one more derivative, it has the 
advantage to involve less terms. Furthermore we do not need to combine with the equation of n. 

Remark 4.3. We note that the profile P is discontinuous at X = 0. Actually this discontinuity is 
compensated by another pressure profile which depends only of t, x. This profile will be constructed 
in the section 14.81 On the other hand we will construct a velocity profile V continuous including 
at X = 0, so that d x P will also be continuous including at X = 0. 

We now use the equation (|132p to get rid of the pressure profile into the equation (|130p . Inverting 
the two sides and dividing by t, we have: 



\ 2 d 2 x V + ^d x V -\v = t(d t V + v° ■ V X V + V ■ V x v° - 2 {V - V f ) - n i 
2 2 nr 



4.3.4 Avoiding a far-field degeneracy 

The vectorfield n may vanish, away the patch boundary, hence so may do the coefficient in front 
of the leading order in the equation (|136p . To remedy to this we consider a function a in the space 

B := L°°([0,T],C°' r (lR d )) nL°°(0,T;C s > r (O±(t))) 
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satisfying the condition infr t T]xR d o =: c > and such that a = \n\ 2 when |<y9°| < tj, and we consider 
for the profile V(t, x, X) the linear partial differential equation LV = where the differential 
operator L is given by L := £ — t(D + A) where £ and A are some operators of respective order 2 
and acting formally on functions V(t, x, X) as follows: 

EV := ad\V + -d x V - -V and AV := V ■ V x v° - 2 {V ' V n. 

2 2 a 

Roughly speaking the equation f)136|) is now hyperbolic in t, x and parabolic in t, X, for t > 0, and 
degenerates into an elliptic equation in X for t = 0. 

The substitution of a instead of |n| 2 is almost harmless since their values are different only for 
l^ ! ^ rj, so that the corresponding values of the eventual solutions V(t, x, ~/=^ ) and V(t, x, ~/=^) 

respectively given by the equations (|136|) and (fTSj) (the equation LV = 0) both tend to as \put 
tends to 0, because of the vanishing condition (I120p for X at infinities. 

4.3.5 Transmission conditions 

Of course we hope to extirpate from the previous equation a non-trivial solution. To this purpose 
an important point is that we look for a solution V with a X derivative dxV discontinuous in 
X = 0. Actually because of the parabolic nature of the Navier Stokes equations, we expect that 
v u and uj u are continuous including through ip° = (these are the Rankine-Hugoniot conditions 
associated to the problem), which lead to the transmission conditions: V and lo° + Cl should be 
continuous, which (taking into account the equalities (|121|) . (|122|) and (|131|) ) is equivalent to the 
transmission conditions: V and n Aw° — |tj| 2 c?x^ should be continuous. More precisely this means 
a priori that 

V\x=o+,<pO =0 + -V\x=o-,cpO =0 - = 0, (137) 
\n\ 2 d x V\ x =o+^o =0 + -\n\ 2 d x V\x= -^o =Q - = -(n A a;°|^o =0 + - n A u\ 0=Q - ). (138) 

Since X is the placeholder for y the function V(t,x,X) needs to be defined only when X and 

v share the same sign. However it is useful to look for a profile V(t, x, X) defined for (t, x, X) 

in the wholedomain T> := (0, T) As a consequence we will actually look at the following 

transmission conditions: for any (t,x) G (0, T) x R : [V] = and [<9v^] = „ — where the 

brackets denote the jump [V] = V\ x =o+ — V\ x =q- across {X = 0} and where are two functions 

in L°° ^(0, T), C s ' r (R d )^ such that u±\o ± ( t } = uj°. More precisely we recall from the compendium 

(precisely the point [6] of Theorem II. lj) that the restrictions of the flow x° of the Euler solution 
on each side of the boundary are analytic with respect to time with values in C s+1 ' r . Thanks to 
Theorem 12.11 there exists some extensions x± analytic on (0, T) with values in C s+1,r (R rf ) of the 
restriction x |o ±- We define the corresponding velocities v± by v±(t, .) := (dtx±)(t, (x±) _1 (^ •)) 
and the corresponding vorticities ui± by ui± := curl v±. As a consequence, with the notation of the 
introduction, u± are in Bp- 

The profile equations (I18p -(|24p announced in the introductory section [l~2l are therefore derived. 

4.4 At initial time 

We expect that the initial values Vo(x,X) := V(0, x,X) satisfy the problem obtained formally 
making t = in the equations (fT8l) - (|24l) . i.e. : 

aod x Vo + —d x V -^V = when ± X > 0, (139) 

l/ | x=0+ - F |x=o- = 0, (140) 
aodxV \ x =o+ - a d x Vo\x=o- = 9o- (141) 
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The functions ao(x) and go(x) which denote respectively the initial value of a and of n A (u^ — tu°_) 
are in C°' r (R d ) n C s > r {0± fi ). 

Proposition 4.1. There exists a unique couple of solutions 

v (x,x) g c°> r (m. d ,p-s(R))nc s > r (o± )0 ,p-s(R)) (142) 

which satisfy the problem U39\ )- (T4U\ )- ljil4i\ )- Moreover 

V Q (x,X) ■ n (x) = for±X>0. (143) 

Proof. Step 1: Reduction. We first reduce the transmission conditions to the homogeneous case 
by defining the functions Vq by 

Vb:=M)±- — {e TX -e^ 2X ) when ± X > 0, (144) 
2 a 

so that the problem (I139p - (ll40p - (ll4ip is turned into the following one (dropping the tilda and the 
index 0) 

a d x V + —d x V - = f when ± X > 0, (145) 

V\ x=0+ -V\ x=0 - = 0, (146) 
d x V\ x=0+ - d x V\ x=0 - = 0, (147) 

with / G C°' r (R d ,H- l (R)) nC°' r (R d ,p- S(R)) nc s > r (o± fi ,H- l (R)) nc s > r (o± fi , P - S(R)). 

Step 2: Cutting-ofF the unbounded coefficient. Because of the unbounded coefficient X 
in the equation ([139p the previous problem does not enter in the classical theory of elliptic problems 
(with x as parameter). To remedy to this we introduce a cut-off. We consider a > and a smooth 
function \<t such that Xa(X) = X for \X\ < a, Xcr(X) = 3a/2 for \X\ > 2a and ||xjr lli°°(R) < 1- 
We will work with the modified equation: 

a d 2 x V + Xa ^ d x V - ~V = f when ± X > 0, (148) 

Step 2: Existence and uniqueness of weak solution. The variational formulation of the 
problem l[I3S]l - ([I3B]l - ljI37j) reads: for / G L 2 (R d ,H~ 1 (R)) find V G L 2 (M d , i7 1 (IR)) such that 
B a (V,W) = - < f,W > for all W G L 2 (R d , if^M)) where < .,. > denotes the bracket of 
duality between L 2 (R d , H~ 1 (R)) and L 2 (R d , H 1 (R)) , and B a (V,W) is the following bilinear form 
on L 2 (R d ,H 1 (R)) x L 2 (R d , H^R)): 

B a (y, W) := / a d x V.d x W - ^^-d x V.W + \v.W. (149) 

JR d xR * * 

Since the bilinear form B a is continuous (thanks to the cut-off x<j) and coercive: 

B a {V, V) = f a \d x V\ 2 + + ^4^)\V\ 2 (150) 

we infer from Lax-Milgram theorem that there exists a unique weak/variational solution of problem 
(fM- (fTi6) -(fT471). 

Step 3: V G L°°(R d , i? 1 (R)). We now turn our attention to regularity restricting ourselves for 
shortness to establish a-priori estimates. First multiplying the equations (|148p by V and integrating 
(only) in X yields for any x G R d 

[ a \d x V\ 2 + \(l + ^l)\V\ 2 dX = - « f,W » (151) 
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where << ., . >> denotes this time the bracket of duality between H 1 (R) and so that we 
get 

\\V\\L°°(R d ,H 1 (R)) < ll/lli°°(]R<i,H-l(]a))- ( 152 ) 

Step 4: V G C s ' r (C± i o, ^(R)). We follow the same strategy than in sectionEjS By Rychkov's 
Theorem 12.11 there exists some extensions a±o in C s,r (R d ) and 

f± £ C s ' r (R d ,H- 1 (R))nC s ' r (R d ,p- 



of the restrictions of ao and / to 0±q. Then we perform estimates of the solutions V± of the 
problems 

a±, d 2 x V± + ^ld x V ± - W± = f± when ± X > 0, (153) 

V±\ x=0+ -V±\ x=0 - = 0, (154) 
d x V±\ x=Q+ -d x V±\ x=0 - = 0, (155) 

for x running in the full range R d through a Fourier analysis. Finally we observe that V\o± o an d 
V± satisfy both the equations (|148j) - (|146j) - (|147|) for x in 0± i q. Proceeding as in step 3 we get that 
they are equal. 

We are therefore reduced to the case where the functions ao and go are respectively in C s,r (R d ) 
and C s ' r (R d ,H- 1 (R))nC s ' r (R d ,p-S(R)). We are going to prove by iteration for -1 ^ I < s that 



ll^ / |lc i - r (M d ,H 1 (R)) ^ \\J\\C l ' r (R d ,H- 

To do this we make use of a spectral localization with respect to x that is we apply the operators 
Aj to the equations (1148|) to get for j > — 1 the equations: 

aodxAjV + ^d x A 3 V - ^AjV = A 3 f + | [a , A 3 ]d 2 x V (156) 
for ±X > and to the interface condition (|146p - (|147p to get at X = 

A i V| x=0+ -A i FU =0 -=0 J , . 

A 3 d x V\ x=0+ - A,d x V\ x=0 - =0. [ > 

We want to show that 

sup^_i 2 j( ' +r) (||A i y|| L oo (M 3 i ^i (R)) < oo. 

To do this we multiply the equation (|156p by AjV we integrate in X over R±, and we sum the two 
resulting equations, noticing that the boundary term produced by the integration by parts of the 
sum of the respective first terms vanish. We thus get for any x £ R d the identity 

I aoldxA^I 2 + ^|A,y| 2 = - « Ajf + \[av,A j \d x V,A 3 V » . 

JR 4 



Now we also have 
so that 

where Ij(x) denote 



/ 



\[ao,Aj]dxV\ x=o+ — 

\[a^A j ]d x V\ 



a \d x A 3 V\ 2 + \\A 3 V\ 2 < ||A,/||^ 1(R) + I 3 (x), (158) 



Ij{x) := [ IdxAjVl.l^A^dxV]. 
Jr 

We will prove the following commutator estimate: 



34 



l](x) 



lf(x) 



Lemma 4.1. For any x 6 M. d there holds 

sup^_ x 2^ l+r h j (x) < WdxVW^^^y \\dxV\\ cl _ ltrnLoo ^ L2{ ^ 

Let us admit for a while Lemma 14,11 and infer from the estimate (|158|) that 

ll^llc i .r(]R< i ,.ff 1 (]R)) ~ WfWc l > r (wi d ,L 2 (R)) + ^ X ^llc i - 1 'TiL°°(R d ,I, 2 (R)) 

so that starting with the case I = —1 -which is a consequence of the estimate (|152|) - the iteration 
can be done till we get V € C s > r (R d , H^R)). 

Proof of Lemma \4-1\ We will consider only j > 0, the case j = corresponding to minor modifi- 
cations of notation and being actually easier. We make use of the paraproduct writing 

[ao, &j\d x V = [T ao , Aj]d x V + (a - T^AjdxV - A,(a - T ao )d x V 
so that Ij = 7- + Ij + 7j where 

7]^) := / |^A^|.|[r aoJ A,-]5xF|, 

= / |9 x A,y|.|(ao-r ao )A^ x F|, 

-- [ IdxAjVUAjiao-Ta^dxVl 

Referring to the definitions (|32p and ()28[) we have 
[T ao) A,-]a x y = ^[S^ia ,A,]A fc 3 x F 

£ 2 3 ^' f (S k ^a (x)-S k ^a (y))h(2^x-y))A k d x V(y)dy. 

Hence 

|[T ao ,A,]9xy| < 2"^ £ ll^fc-iaolUip / 9 (2^(x-i/))|A^x7(y)|2 3 ^y 
fc>l,|fc-j|<5 ,/rS 

where g denotes the function g(x) := which is in L 1 (M 3 ). Using Fubini theorem we get 

that for any i£R 3 

I)(x)< Y, 2 ~ J I 9(2 j (x-y))( f \d x A J V(x,X)\.\d x A k V(y,X)\dX)2^dy. 

We now use Cauchy-Schwarz inequality to get for any i£R 3 

Ij(x) < Y 2^\\d x A j V(x,.)\\ L , m [g(2^x-y))\\d X A k V(y,.)\\ LHR) 2^dy 
fc^l,|fc-j|<B JU3 

< Y 2~ : >\\d x A j V(x, .)\\ L 2 {R) \\d x A k V(y, .)IIl-(r3 iL 2 (]R)) 

fc>l,|fe-j|<5 

hence 

BUp^ 2 2 ^lj(x) < PxV^^^yPxV^^^ 

Similar bounds for 7? and 7? can be obtained, mixing once again classical paradifferential arguments 
with a Fubini argument. 

□ 
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Step 4: V G C s > r (R d ,p - <S(R)). Mimicking Step 3 we can prove that 



Now we are going to prove by induction on I and then by induction on k that for all k and I in N 

V kjl>0 := X k d l x V € C s > r (R 3 lP - H X (R)). (160) 

Thanks to Sobolev embeddings we will infer (|142p . The point is that Vk,u let us simply write V, 
verify a problem of the form 

a d 2 x V+ ^d x V- ~V= F when ± X > 0, (161) 

V\x=o + ~ V\ x=0 - = J, (162) 
aodxV\ x =o+ ~ a dxV\ x =o- = J, (163) 
V{x,X)^0 when X -> ±oo. (164) 

where JP (respectively J and J) is in C s ' r (M. d ,p — L 2 (M)) (resp. in C s ' r (lR a! )) thanks to the previous 
steps. Once again we can reduce the transmission conditions to the homogeneous case by defining 
the functions V by 

V :=V±-((2J + -)e^ x -(J + -)e^ 2X ) when ± X > 0. (165) 
2 a a 

which satisfy a problem of the form 

a d 2 x V + jd x V - W = F, (166) 

V\ x=0+ -V„\x=o- =0, (167) 
ao<9x^|x=o+ - ao9xV| x=0 - = 0, (168) 
V(x,X)^0 when X -> ±oo. (169) 

where F is in C s > r (R d , L 2 (R)) so that we can apply (fT59l) . 

Step 5: y G C°' r (IR d ,p — 5(R)). It is actually given by the analysis of step 4. 

Step 6: Passing to the limit in a. We finally let a goes to infinity. The estimate (|177p is 
uniform with respect to a. Using weak compactness we can pass to the weak limit so that we get 
a solution to the uncut equation. 

Step 7: proof of (11431) , It suffices to use uniqueness argument based on equality (I15ip . 

□ 

4.5 Existence and uniqueness of the layer profile 

Let us now study the time-dependent equation (|18p . We denote by T> the domain T> := (0, T) 
R, by T>± the restrictions T>± := (0, T)xR d x R± on each side of the boundary T := (0, T) x R d x {0} 
and by A, L and L the operators of respective order 0, 1 and 2 acting formally on functions V(t, x, X) 
as follows: 

AV: = V.V^-^ V - V ^ n n 

a 

LV := jd x V ~\v- t(d t V + v° ■ V X V + AV), 

LV := ad x V + LV. 
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The profile problem now reads as follows 

LV = f on V ± , ([V], [d x V]) = (0, 5 ) on T, (170) 

the decreasing of V for large X being encoded in the choice of the space E\ := L 2 ((0,T) x 
R d , H 1 ' 



Theorem 4.1. For any f G E' 1; for any g G L 2 ((0,T) x R rf ) i/iere exists exactly one solution 
V £ E\ of \17(ty . In addition the function y/i\\V(t, ., -)llL 2 (R d xR) * s continuous on (0, T). 

The first equation of fll70|) is satisfied in the sense of distributions on both T>± and the sense 
given to the jump conditions is explained in the proof below. 

Proof. We consider a > and a smooth function Xa such that Xa{X) = X for \X\ < a, x<r(X) = 
3a/ 2 for \X\ > 2a and ||XctIU°°(r) < !• We will work with the modified operators L CT := L + 

x ° <yX ^~ X dxV and := ad\V + L a V whose coefficients are bounded. Even more the coefficients 
of the first order part are Lipschitz. 

We now explain the meaning of the jump conditions in the equation (|170p . Since V is in E\ 
the jump [V]|r is in L 2 (T). To give a sense to the jump of the derivative [<9x^] we will use the 
equation. For any V in the space E 2 := {V £ Cq{V)/ V\t> ± G C°°} and W in H 1 ^) we have, 
integrating by parts, the following Green identity: 

V / L a V.W = - [ ad x V.d x W+ [ L* a W.V - [ a[d x V].W-T f w.V (171) 
± Jv± Jv Jv Jr Jf 

where T := {T} x R d x R and where L* denotes the operator (the adjoint of L a ) 

KV ■= -^4p-dxV - i(l + xUX))V + t(d t V + v° ■ V X V + (1 + divv - A)V). (172) 

In fact less smoothness is needed. Let us introduce the Hilbert space E4 := {V G E\j L a V G 
i7 _1 (D)} endowed with the norm [|V||.e 4 := HVH^ + II-^ct^IIh- 1 ^)- Thanks to a classical lemma 
by Friedrichs |52j the space E2 is dense in E4. 



Lemma 4.2. The map 



VeE 2 ~t:={ " ld T l V] ° n f r (ITS) 



extends uniquely to a continuous linear map from E4 to H 2 (Y U f ) and Green's identity {111 ) is 
still valid for any couple (V,W) in E4 x H l (T>) in the generalized sense that 

< L a V,W > H - lmm{v) = jjVllW - ad x V.d x W)- < r,W\ ruf >^ (Rjf) ^ (Rjf) . (174) 

Proof. Let V be in E2 and W be in H^(T U f ). There exists a function W in H l (T>) such that 
W^lruf = From Green's identity (|17ip we infer that 

I | ru _ rW\ < C\\V\\ E4 \\W\\ HHV) < C||^lbJ|W||^ (ruf) . 

Hence by Hahn-Banach theorem we get the existence of a continuous extension, which is unique 
because of the density stated above. □ 

We therefore have given a meaning to the problem (|170|) , This meaning can seem weak but the 
next result said that is actually strong. We introduce the space H 1,2 (T>) of the functions V G H l (T>) 
such that d x V\x> ± are in L? . 
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Lemma 4.3. If V £ E4 satisfies the jump conditions [V] = and [dxV] = g on T in the sense 
given by Lemma \4-%\ then there exists a sequence V £ in H 1,2 (T>) converging to V in E4 and a 
sequence g e converging to g in L 2 ((0, T) x M. d ) such that [V £ ] = and [dxV 6 ] = g £ on V. 



Proof. As this kind of process is very classical, see for instance Rauch [99], we only briefly sketch 
the proof. The idea is to construct the sequence V s by convoluting in the variables t, x only to 
preserve the jump conditions, to use Friedrichs lemma to prove the convergence in E4 and then to 
gain the extra X derivative, that is to prove that the V s are in H 1,2 (T>), thanks to the equation. □ 

We will now prove uniqueness as a consequence of the following estimate: for any function V 
in Ex satisfying 

L a V = f when on V±, [V] = and [d x V] = g on T, (175) 

there holds 

\\V\\ei ^ ll/lls; + IMIl^xr)- (176) 

Let us first see the case where V is in H 1,2 (T>). Then using Green's identity (I17ip with V = W 
and density we get that 

/ \V\ 2 + \d x V\ 2 + T [\V\ 2 < [ t\V\ 2 + \f.V\ + / \g.V\, 
Jv Jf Jv Jr 

actually not only for T but for any t in (0, T), so that using a Gronwall lemma we get 
\\V\\ El +y/f\\V(T,.,.)\\ 

L 2 (R d xR) < \\f\\E[ + IMlL 2 (R d xR)- ( 177 ) 

Now to deal with the general case we introduce the Hilbert space E3 := {V 6 E\j L a V 6 E[} 
endowed with the norm [|V||.E a := + [| Argils' • From its definition (I172() we see that L* a V is 

in E[ whenever V G E-$. Estimate (|176p will be deduced from: 



Lemma 4.4. The map 

(V,W) £E 2 x E 2 ^ p:= 

extends uniquely to a continuous bilinear map from E3 x E3 to (Lip(T U T)' and Green's identity 



a[d x V].W on T 
TV.W on f 



7771) is still valid for any couple V, W in E3 x E% in the generalized sense that 



< L a V, W > E ^ E = - jf ad x V.d x W+ < L%W, V > E ^ El - < p, 1 > (Lip(ruf) ^ K ruf) ■ (178) 

Proof. Now let V, W be in E% and let ip be in Lip(r U T) there exists a function in Lip(2?) such 
that 3>| ru f = <p and such that ||$||Lip(x>) ^ IMlLip(ruf) w ^h c independent of ip. From Green's 
identity (jl71l) we infer that 

f p<p = - f L a V.(pW- [ ad x V.d X (pW+ [ L* a (pW.V 
Jr Jv Jv Jv 

so that 

\J M < \\nu P {vy\\LaV\\ E[ -\\W\\ El + ||$|| L ip(2>)-l|y|k-l|Wlk + \\LlipW\\ E[ .\\V\\ El {l^) 

< CMk^vyWWWW (180) 

Hence by Hahn-Banach theorem we get the existence of a continuous extension, which is unique 
because of the density of E 2 in £3. □ 
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In order to prove the existence part of Theorem 14.11 we shall need another Green formula which 
involves the complete transposition of the operator L. At a smooth level that is for any V, W in 
the space E2 this Green identity reads: 

f L a V.W = I L*W.V- i a[d x V].W + f aV.[d x W] - T I [w.V (181) 

where L* = ad\ + L* denotes the adjoint of the operator L a . 
Lemma 4.5. The map 

V.E^r^l^Y ( 18 2) 

extends uniquely to a continuous linear map from E4 to H~ j(ruf) and Green's identity M81\) is 
still valid for any couple V,W in E4 x H > 2 (T>) in the generalized sense that 

< L a V,W > H - HV) , m(p) =< L:w,V> E ,, El - < r,W\ ruf > H -^ ruf)tHkrut) (183) 

+ f[d x W).aV\ T . 

Proof. By adapting classical lifting method we get that for any W\ in H?(T U f) and W2 be in 
L 2 (r) there exists a function W in H X > 2 (V) such that W\ Vvjt = Wi, [d x W]\ r = W 2 . The proof 
then follows the same lines than the one of Lemma 14, 4( we use Green's identity (|18ip for V be in 
E2 and a W obtained by lifting W\ in H 2 (T U V) and W2 = 0, next we use Hahn-Banach theorem 
and density. □ 

Proceeding as previously we get that the estimate (|176p also holds for the adjoint operator. We 
infer the existence of a solution to the direct problem by using Riesz's theorem and Green's identity 
dl83|). 

We now use again a sequence of approximation of the solution, and use the linearity of the 
problem together with the a-priori estimate (|177l) to show that this sequence is a Cauchy sequence 
in the space of the functions V such that V^II^C^ •> •)llL 2 (R ti xR) i s continuous on (0, T). Completeness 
yields the conclusion. 

We finally let a goes to infinity. The estimate (|177p is uniform with respect to a. Using weak 
compactness we can pass to the weak limit so that we get a solution to the uncut equation. 

□ 



4.6 Smoothness of the layer profile 

We are now studying the smoothness of V(t, x, X). It is convenient to recall here a few notations 
given in the introduction. For any Frechet space E of functions depending on t, x and possibly on 
X we will denote Ed the set 

Ed := {/ eE/3C> 0/ (§]^)keN is bounded in E}, 

where D denotes the material derivative D := dt + v° ■ V x . The spaces denote 

A := L oo ((0,T),C ' r (lR d ,p-5(lR ± ))) n L°°((0,r), C a - r (0±(t),p - 5(M±))) , 
B := L°°((0,T),C ' r (lR d )) n L°° ((0, T), <T' r (0 ± (t)) . 

Moreover we recall from section 14.3.51 that the functions a and A are in Bd- 
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Theorem 4.2. For any f £ Ad, for any g\,g 2 6 Bd there exists one solution V £ Ad of 

LV = f on V±, ([V], [d x V]) = (g u g 2 ) on V. (184) 

Proof. We first reduce the transmission conditions on the internal boundary T to homogeneous 
ones by defining the functions V by 

V :=V±^((2g 1 +g 2 )e fX -( gi +g 2 )e^ 2X ) when ± X > 0. (185) 

so that the problem now reads as follows 

LV = f on V±, [V] = [8 X V] = on T, (186) 

with / 6 Ad- Let us stress that there is no loss of regularity in this lifting process, since the 
X-derivative has been applied innocuously to the second term in (|185|) . as well as D since gi,g 2 
and a are in Bd 

We first prove that V is in C D where we denote C := L°°((0,T) x R d , iT^R)). Let us first 
establish an a-priori estimate, applying for any fc£N the field D k to the problem (I186p to get 

(S - k)VW = /M + + tAVW on V ± , [V^] = [d x V^] = on T, (187) 

where we denote by := D k V the kth iterated derivative of V along D and where := 
Ysi=i fi i where /{ := D k f whereas and fjj^ denote respectively the commutators: 

fT:=[D\S]=J2 (^)D k - l a.dlV^\ ff ] := -[D k , tA] = - g QlJ^tA.^. 

The last sums have to be omitted when k = 0. Here we have used that [D k ,tD] = kD k . We now 
multiply the first equation of (|187|) by U^' and we now integrate w.r.t. X only. This yields for any 
t,xe (0, T) x M d the estimate: 

i\d x V [k] \ 2 + {k + -) [ \V [k] \ 2 ^\ f / t|V-[* +1 I.V-[ fc I|+Ci / t\V [k] \ 2 . (188) 

4 J_r Jr Jr Jr 

Using the condition (|23p the l.h.s. of (|188|) is larger than 

of \d x v^\ 2 + -(k + i) [ \v^\ 2 . 

Jr 4 J R 

We now bound the r.h.s. of (|188p . Using Cauchy-Schwarz and Young inequalities we have 

f i/w.ywi ^ 4 i |/[fe]|2 + k + i r ! - W|2 _ (189) 

JR « + 1 JR 4 JR 

Integrating by parts yields 

k— 1 

I / jf-t^l < £ (f) / l^^9xU H 9xU W |. (190) 

JR ^ \V JR 

Since a is analytic, there exists C a > s.t. for any Z G N ||D'a||s ^ C^Z! so that using Cauchy- 
Schwarz and Young inequalities yield 

\!ft ] .vw\ < c^E ^ci-'na^ w ll) 2 + | /W w l a > (i9i) 

JR l=Q L - z JR 



40 



where we denote here ||/|| := (f R \f(t, x, X)\ 2 dX) 2 . In a similar way there exists C%,Ca > s.t. 

\[ fl k] .vW\ < C 3 C£ jC k A ~ l \m\) 2 + l [\V [k] \ 2 . (192) 



z=o 

Finally for < t < T we have 



^ k±l f |v-[*]|2 + _^_ I |y[fc+l]|2_ (lg3) 

8 /m A; + 1 



Hence 

fe-i 



A;— 1 

+ c 3 (y ^c k A - l \\v^\\) 2 + [ \V&+Y + CiT[ \vw\ 



=0 

thus we infer -keeping the notations C4-C3 for their squareroots- that 

|||^VW|| + ^||VW||< jj=\\fW\\ + C 2 J2 jCt l \\d x V^\\ (194) 



1=0 

We introduce the functions 



k— 1 / 



||F [fe] || \\d x V [k] \, 
a k (t,x) := —t , 6fc(t,x) := — — —77== and fk(t,x) 



k\C k ' fc!C fe Vfc + l (A; + l)!C fc ' 

where C is a positive real which will be chosen in a few lines. Dividing the estimate (|194h by 
k\C k y/k+l yields 

k—1 k—1 

£& fc + I Qfe *:2/ fe + C 2 ^ (^) fc A + C 3 ^ (^) fc -% + 74rCa fc+1 + C lv / Ta fc . (195) 

1=0 1=0 

We choose C large enough so that max( c 2 , J^ 3 ) ^ min(|, A) and then T is chosen small 

enough so that ^f^TC ^ ^ and that C\y/T_ ^ gj. Hence summing over k S N the estimates (I195D 
yield the a-priori estimate: for any t, x € (0, T) x R d 

fceN feeN 
We now define the iterative scheme (V n ) n ^ by setting V° as the solution of 

£V° = f on X>±, = = on T, 

and V n+1 as the solution of 

£ yn+i = f + t ( D + Qn p±j [yn+ij = [a x y n+1 ] = on r. 
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A moments thought convinced that V is in Co- Then proceeding as in the proof of the estimate 
(|196p we infer the convergence of the iterative scheme for any t, x € (0,T) x R rf to a solution V of 
the problem (|186H . Using several time slices yields that V is in Cd- 

Now to prove Theorem 14.21 we increase the smoothness with respect to x thanks to the oper- 
ators Aj of spectral localization. We proceed as previously dealing with the additional spectral 
commutators as in section |4~41 (in particular using the commutator estimate (I4.ip ). 

If we denote by B the space of the functions f(t, x, X) with B smoothness in t, x with values in 
iJ x (R), this yields that V is in Bf). Then we prove by induction that X k V is in Bd for all k in N. 
Finally we use the equation to increase by induction the number of derivatives with respect to X 
and get that V is in Ad- 

□ 



4.7 Other properties of the profile 

Let us now prove that the normal projection V ■ n vanishes. We multiply the equation (|130p by n 
and take into account the equation (jHj) for n(t, x) so that we get 

\n\ 2 d 2 x (V ■ n) + ^d x (V ■ n) - ~(V • n) = t(d t + v° ■ V X )(V ■ n). (197) 

Moreover taking also the scalar product of the transmission conditions on the internal boundary 
X = with n we get 

(V-n)\ x=0+ = (V-n)\ x=0 -, (198) 
d x (V-n)\ x=0+ = d x (V-n)\ x=0 -. (199) 

Proceeding as in section 14.51 we get an energy estimate for the problem (|197j )- (jl98p - (|199p from 
which we infer that V ■ n vanishes. 

When s > 2 the first and second time derivatives of the functions n and ojq are in L°° ( [0, T] , C 0,r (M. d ) ) . 
This allows to get some estimates for the second time derivative d 2 V by commuting d 2 with the 
equation (|18p - (|24| ). Hence we can define a trace at t = of all the terms in the equations (|18|) -(|24| ) . 
We hence get that the trace V\ t=0 satisfy (fl^ - (fT4T)]) - (fT4T]l . 



4.8 Higher order expansions 

We are now concerned with the following terms in the expansion with respect to vt of the solutions 
of the Navier-Stokes equations. Actually in this section we show that if the initial data is piecewise 
smooth on each side of the interface {<^° = 0} -that is if s = +oo- then it is possible to write a 
complete asymptotic expansion of the vorticity of the form: 

u v (t,x) = YV^ j SP(t,x, ip0 ^ )+O{Vri° o ), (200) 

where the first profile Q° is the one constructed in the previous section: Q° := Q. 

Let us explain the underlying intuition to guess the expansion (|200p : when the initial data 
is globally smooth, the vorticity tu u given by the NS equations admits a regular (Taylor-like) 
expansion: 

w"(t,x) :=J2(vt) j Uj{t,x) + 0({vt)°°), (201) 

so that we expect that when the initial data is piecewise smooth on each side of the interface 
{(f° = 0}, the same should still hold true when we incorporate the fast scale - so that we 
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should write uj u through the ^-dependent profile: 



LO v (t,x) =n u (t,x, (202) 



where Q u (t,x,X) admits a regular expansion: 

n u (t, x, X) := ^ x, X ) + 0(Vrt°°). (203) 

Once again to prove this we will consider the velocity formulation, looking first for a determination 
of the velocity profiles of the expansion 

v"(t,x) = v°{t,x) + V vV V j (t,x, 99 °^ ) ) + 0(v / ^° O ), (204) 

where the profile V 1 is the one constructed in the previous section: V 1 := V. Then we will recover 
the vorticity profiles by plugging (|200p and (|204p into the relation (fl~3"|) and equalling the terms of 
order \fvt J '; what leads to: 

= curUW + n A 9 X V J ' +1 . (205) 

To guess the amplitudes in the expansion of the pressure we have to remember that we follow a 
Rankine-Hugoniot(-BKW) approach looking for profiles such that the expansions solve -formally- 
the equations (JI])-© on each side {±y?° > 0} plus the continuity of the velocity, of the pressure and 
of the vorticity on the interface {ip° = 0}. Taking the divergence and the normal scalar product 
yields the pressure problem: 

Ap v = -div (v u ■ Vv u ) (206) 
bl = 0, (207) 
[8 n tf\ = [(_„" • Vv u + pAv").n], (208) 

where the notation [•] stands for the jump across the interface {ip° = 0} (that is for a function 
f(t,x): we denote [/] := /Lo = o+ — /Lo =0 -)- If we plug the expansion (|204p into the r.h.s. of the 
equations (I206p - (|208p we are led to look for a pressure expansion of the form 



jfft x) = p°(t, x) + Y J ^~ P J (t, x, + 0{^T), (209) 

The profiles above are of the following form: for ±X > 0, 

U(t,x,X) := U(t,x) + U(t,x,X), (210) 

where the function U(t,x,X) is rapidly decreasing when ±X — ► oo, and the letter U is the place- 
holder for the IP , the V 3 and the P 3 . We will refer to the term U as the regular part and to the 
term U as the layer part. The layer part P 2 is equal to P 2 = tP where P is the profile of the 
previous section. This possibility to be smoothly factorized by t is very particular to the order 2. 
It comes from the fact that the (a-priori) main contribution given by the Laplacian in the r.h.s. 
of the equation (j208H vanishes thanks to the orthogonality property (11311) . It is therefore fair to 
acknowledge that the heuristic argument given in section 14.21 to guess the first order expansion 
01241) is at less not very robust and leads to a correct amplitude almost by chance ! Furthermore 
to satisfy the pressure continuity we will have to add to this layer part P 2 a regular part P 2 , as 
anticipated in Remark 14.31 
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To determine the velocity and pressure profiles we proceed iteratively, determining at the step j 
the velocity profile V 3 , the regular part of the pressure profile P 3 and the layer part of the following 
pressure profile P 3+ , from the profiles already known by the previous steps. The step j = 1 was 
done in the previous section. We now explain how to do a step j ^ 2 when the previous ones are 
done. 

Remark 4.4. Regarding the transmission aspect of our strategy we refer to [105J, in the setting of 
the approximation of semi-linear symmetric hyperbolic systems of PDEs by the vanishing viscosity 
method. We also mention [70], [7_l] by Gues and Rauch in the context of internal waves for semi- 
linear symmetric hyperbolic systems and [73 J by Gues and Williams in the context of viscous shocks 
profiles. 

4.8.1 Jump conditions 

The transmission conditions read: 

V J \ x =o+,<p°=o+ = V 3 \ x =o- ;tp o=o-, (211) 
(ad x V 3 +d n V 3 ' 1 )\ x=0+tV o =Q+ = (ad x Vi + d n Vi- 1 )\ x=0 - yip o =0 -, (212) 

P J \x=0+,tp°=0+ = -P J |x=0-,y?°=0- (213) 

We use the vectorfields w k defined by (|4"Tj) to decompose the vectorfields as u = ^(u.n)n+u tan where 
Utaxi = Ck(V.w k )w k . In what follows we will use the notation [•] for two different meanings: 

[U] ■= EU=o+ - ZZV =0 - and [U] := U\ x=0 +^o =0 + - U\ x=0+>(p0=0 +. (214) 

By assumption the previous velocity profile V 3 satisfies 

V 3 ~ 1 \ x =o+ ,<p°=Q+ = V j ~ l \x=o-,<pO=o-- (215) 
As a consequence there holds 

[dn^ti; 1 ] = [(^' _1 )t«J, and [{d n V 3 ~ l ).n]n = [^((V^-^nJn.)], 
and the jump conditions A211 |> - (|212[) - fl213[> reduce to 

[V 3 .n] = -[V j .n], (216) 

[*£nl = "Kan], (217) 

WxVH = ~[dnVL\ (218) 

[adx{V j .n)\ = -[d n {V 3 -\n)l (219) 

[P 3 ] = -[pi]. (220) 

4.8.2 Divergence equation 

First plugging the ansatz (|204p into the equation ([2]), equalling the terms at order \J~vt 1 yields: 

d x (V j .n) = -dw x V 3 '- 1 . (221) 

We first begin to take the limit ±X — > oo of the equation to get 

diVa-V? -1 = 0. (222) 

Then subtracting the equation (|222p to the equation (|22ip yields 

d x (V j .n) = -dw x V 3 '- 1 . (223) 
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For ±X > 0, we integrate between X and ±00 to find 

/•±oo 

V 3 .n :=- div x V 3 ' 1 . (224) 
Jx 

This already yields that ()219j) is satisfied. To do this we will use that the divergence can be rewrited 
as 

d\v x V = -d n (y.n) -J2 a kW k -Vx(w k .V) 

k 

where the w k are tangential to the hypersurface and the a/% are scalar. Since &\v x V_ j ~ l = there 
holds 

[d n (V j -\n)} = [a^a k w k .V x (V j -\w k )} 

k 

= -[a^a^.V,^" 1 ./)] 

k 

because of (|21ip . Then 

k 

= [adiv,^'- 1 ] 
= —[adxVKn] 

because of (12231) . 

Moreover proceeding in the same way at the next order yields 

dW x V j = 0. (225) 

4.8.3 Velocity equation 

When plugging the ansatz (|2U4p and (|209|) into the equation ([T]) the term 

that should a priori appears in the equation, vanishes thanks to the eikonal equation satisfied by 
ip°, and equalling the terms of order ^pvt 1 one get one equation of the form 

(L - J -^~)V j = P + V x P j + d x Pi +1 n, (226) 

where f 3 would be decomposed into f 3 = f J + f 3 -as the profiles U were in (I210p - and contains 
terms already determined so that it could be considered as a source term at this step. 

4.8.4 Regular part 

We first begin to take the limit ±X — > 00 of the equation (I226|) in order to kill the layer parts and 
to obtain the following equation for the regular part: 

tv 3 + t(d t V j + v° ■ V X V> + V j ■ V x v°) + V X P> = I 3 , (227) 

Since the hypersurface {t = 0} is characteristic the equation (|227p does not need any initial condi- 
tion. 
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Lemma 4.6. There is only one couple (V_ 3 ,P 3 ) piecewise smooth satisfying the equations (227)- 
1225\) on each side 0±{t) of the interface, and the transmission conditions 1216\) - [22U\) . 



Proof. The idea is to recast the problem into an equation of the form 

|v? + (dt + v° ■ V x )tV j = IICV?) + fj- (228) 
on each side 0±(t), where II is a pseudodifferential operator of order 0. For instance if we denote 

g ■= - 3 -¥i - t(d t v j + v° ■ v x Vi + y? ■ v x v°) + p. 

Taking the divergence and the scalar product by n of the equation (|227f) we get the following 
problem for the pressure: 

AP j = dive/ for x £ 0±(t), (229) 

[d n P j }= = [g-n], (230) 

[pi] = -[pi], (231) 

Using that the regular parts of the velocity profiles -that is here v°, V? and Vj 7-2 - are divergence 
free (see the conditions ([2]) and (|225p ) we simplify the equation (|229[) into 

= div(f -V3-V x v°)-^2d k vldiVi. (232) 

i,k 

To simplify the equation (|230p we apply the operator dt + v°- V x to the condition (|216f) and simplify 
thanks to the equation of n. This leads to 

[d n P 3 ] = + t{(d t + v° ■ V x )Vi + V? ■ V x v°) + f].n (233) 

These simplifications done it becomes clear that the operator Ii{V^) := —V X P? — V_ 3 • V x v° is of 
order 0. 

We now set 7 := |— 1 and we denote 0± the space-time domains 0± := {(t, x) £ (0, T) xW 1 / x £ 
0±(t)} and H the operator: 

H := 7 + (d t + v° ■ V x )t. (234) 

We have 

Lemma 4.7. Both problems 

HV = l for Q± (235) 

are well-posed: for any f in L 2 (0±) there exists only one solution V_ in L 2 {0±). Moreover V_ in 
C((0,T),L 2 (O±)). Iff is in L°°({0,T),C s ' r (O±(t))) then V in L°°{(0,T),C s ' r {O±(t))). 

Proof. Existence and uniqueness in L 2 are based as in section 14.51 on a Green identity: for any 
smooth function W_ there holds 

f HV.W = j [ V.W- f tV.(d t + v° -V X )W + T [ V.W. (236) 
Here we have used that 

d t f = [ (d t + v°-V x ). (237) 
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This identity leads to the a-priori estimate: 

(7 + h [ V 2 + \ I Y- 2 < I f- W 
1 Jo± z Jo±(t) Jo± 

This estimate, by a weak-strong argument, leads to uniqueness. Besides by Riesz theorem we get 
the existence part. Using Rychkov's extension theorem 12.11 the estimates reduce to estimate in the 
full plane for globally smooth coefficients and source. Then we apply the operators Aj to get 

HA j V = l + t[v°-V x ,A J ]V. (239) 

We multiply this equation by p\AjV\ p ~ 2 .AjV , where p > 2, and we integrate by parts in t,x over 
(0, T) xM. d . The commutators are estimated through paraproduct estimates which follow the same 
lines as the one of the previous sections. Finally we let p goes to infinity. □ 

From the previous lemma we infer the existence and uniqueness in I? for the equation (12281) . 
To get piecewise smoothness we notice that the operator II satisfies the transmission property so 
that it preserves piecewise smoothness at any order (we recall that in this section the interface is 
assumed to be smooth so that we only need here to refer to classical works as for instance [100J ) . □ 

4.8.5 Layer part 

The tangential part of Y3 is i n general discontinuous at {<^° = 0}. This jump will have to be 
compensated by the layer part V 3 , which we now look for. We subtract the equation (|227p to the 
equation (I226P to get for V 3 an equation of the form 

(L - 3j ^~)V ] = d x P j+1 -n + f. (240) 

Lemma 4.8. There exists a couple (V 3 , P 3+ ) piecewise smooth satisfying the equations {240\ )- 
{221$ on each side ±X > of the interface, and the transmission conditions {217\ )- (2T8\) . 

Proof. Since the equation (|224p already determines the normal part V 3 .n of the velocity profile V 3 
we will reduce the problem into a problem for the tangential part Vj^n- The pressure P 3+1 can be 
determined in function of V 3 : 

\n\ 2 d x P j+1 = (L- J -^-)(V j -n) + 2tV j ■ V x v°.n - f.n. (241) 

Here we have taken the scalar product of the equation (I21ip with n and use the equation of n. 
Taking (|24ip into account the equation (|240p now reads 

(L ~ J ~=^)V j = "((£- J -^)(V 3 .n))n + 2t(V 3 ■ V x v°.n)n + # an . (242) 
But using again the equation of n we have 

(i-^yVtln = (L-^)V 3 -±((L-i^)(V 3 .n))n 

+ t -{V 3 .n)(-(n.V ' x v° .n)n - l V x v° .n)) (243) 

Mixing (pi2l) and (f2i3l) leads to 

{L- 3 -^){vL) = f %an + tAV 3 (244) 

where AV 3 denotes the local 0-order operator 

AV 3 := 2(V j ■ V x v°.n)n + -(V j .n)(-(n.V x v° .n)n - *V a w°.n)). (245) 

a a 

Thanks to the analysis of section [431 there exists V^in satisfying the equations (|244p - (|217p - (|218p , 

□ 
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4.8.6 On the smoothing of the level function 

We should be tempted to take into account the viscous smoothing of the level function ip° that is 
instead of the solution tp° of 

Dv° := d t if° + v° ■ V x ip° = 0. 

to consider the solution (p u of the "viscous eikonal equation": 

d t ip v + v v • Vsp" = vA x tp\ (246) 

with, again, the initial data: ip u \t=o = (fo- Actually reading again the formula (|126|) - (|127|) - (|128|) we 
notice that these terms are the ones in prefactor of -j^dxU when one applies the transport-diffusion 

operator to a function of the form U(t,x, ^=^). 

In that case proceeding as previously we have that (p v admits an expansion of the form 

^(t,x) ~ rt(t,x) := <p%x) + irt*(t,x,^f-) (247) 



where he profile &(t,x,X) is defined as §(t,x,X) := &(t,x) + &(t,x,X), where the $ are the 
solutions of 

(1 + tD)<t = A x cp° (248) 

and $ is the solution of the equations 

(S - - - tD)® = 0. (249) 
with the following transmission conditions across X = 0: 

[*] = - t-), and [d x $] = 0, (250) 

where $_|_ are suitable extensions of Let us show that A x Lp v a does not have a jump. We have, at 
least in a neighborhood of {<£> }, the following identity: 

A x .^ = A x <p° + ad 2 x § + Vvt{A x <p°.d x $ + In ■ V x d x $) + vtA x §. (251) 
In the expression above the profiles are evaluated in X = y . Using the equation (|249|) we have 

[ad x $\ = [(1 + iD)*]. (252) 

Since D := dt + v • V ' x is tangent to the hypersurface {<p°(t, x) = 0} we infer from (|250p that 

[ad x $] = -[(1 + tD)$\. (253) 

Using the equation (|248|) we get [a9^$] = — [A x tp°], so that 

[A x <p» a ] = VVt[(A x <p°.d x $ + 2n • V x d x *)] + vt[A x $]. (254) 

Besides differenting with respect to x the transmission conditions (|250p we get [A^^] = 0. 

One can see the link between the two points of view by a Taylor expansion: for any smooth 
profile U there exists a smooth profile U' 9 such that 

U(t, x, = U(t, x , + WtU\t, x, (255) 

yjvt Jut Jvt 
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In effect it is sufficient to define U b by U b (t,x,X) := U*(t,x,X,$(t,x,X)) where U^{t,x,X,h) := 
h- l (U\t,x,X + h)-U\t,x,X)). 

We believe that the point of view of the phase smoothing adopted in this section could be 
interesting in view of future extensions to singular vortex patches. Let say for instance that V^o 
has a jump on the hypersurface tpQ. Then the corresponding solutions <p v given by (j246H admit an 
expansion of the form: 

<f?(t,x) ~ </P(t,x) + Vuiip 1 (t,x, - ^1^ )) (256) 



where tp u also solve the eikonal eiquation: 

W + v° ■ V X *P P = i/A^, (257) 

^"|t=o = ^o. (258) 

Let assume that ipo is smooth so that ip u is well-approximated by the solution ip° of the inviscid 
eikonal equation: 

d t ip° + v° • V X V° = 0, (259) 

V>°|t=o = ifo. (260) 

As a consequence the expansion (|256h can be simplified into 

ip v (t, x) ~ p°(t, x) + vW(*, x, ^-M> )). (261) 



Now we plug this later expansion (I26ip into the expansion (I262p to get 

^ft*)^ (262) 

where 

n(t,x,Xi,X 2 ) :=0(t,x,Xi +< / 9 1 (t,x,X 2 )) (263) 



4.9 Tangent parallel layers do not interact 

As mentioned in Remark 11.31 the proof a la Chemin of Theorem 11.11 also succeed to cover the case 
where the initial vorticity ujq is discontinuous across two hypersurfaces {ipo = 0} and {(fo = rj}, 
where r/ > 0: it yields that the corresponding solution of the Euler equations has -at time t- 
a vorticity ujq piecewise smooth and discontinuous across the two hypersurfaces {cp°(t,.) = 0} 
and {(p°(t,.) = 77}, where ip° is again the solution of the transport equation ©-((HI). For the 
corresponding solution of the Navier-Stokes equations two layers of width \pui develop around the 
hypersurfaces {(p°(t, •) = 0} and {<p°(t, .) = ry}. When t proceeds it comes a time when the layers 
overlap. However they do not interact and the NS velocities can be described by the superposition 
of the two layers that is by an expansion of the form 

v»(t,x) ~ v°(t, x) + VUt( V(t, x, + W(t, x, V ^- V )) , (264) 

V yvt \Jvt > 

where the profile V is again the solution of the equation (fT8]) with the transmission conditions (j2"4"j) , 
whereas the extra profile W is the solution of the equation (|18|) . with the transmission conditions: 

[W] = and [d x W] = - nA( "°-" 0) , 

a 



49 



where the brackets denote the jump of W(t, x, X) across {X = 0} that is [W] := W|x=o+ — VF| v=o- 
and tv± are some well-chosen extensions of the restriction of the Euler vorticity to both sides 

{:(/('••; '/) >o}. 

The point is that both profiles V and W satisfy the orthogonality condition V ■ n = W ■ n = 
so that the Burgers term 

{V + W) -nd x {V + W) 

which should induce a nonlinear coupling of the two layers identically vanishes (as the self-interaction 
did in section I4.3.2[) . 

In the same vein -using locally this argument- we can infer that there is no non-linear interaction 
between two vortex patches tangent in one point. 

4.10 Well-prepared expansions 

Let us also stress that it is also possible to construct some asymptotic expansions of the form: 

cj»(t,x)=J2VJ W(t,x,^^-) + 0(Vt°°), (265) 

v u (t, x) = v°(t, x) + Y] V j (t, x, fJhA) + 0{y/v°), (266) 
p»(t, x) = p°(t, x) + Y, >/? P J (t, x, + 0(V^°), (267) 

for both the Euler and the NS equations. Let us see how the construction is modified, at least for 
the first velocity profile. 

For the Euler equation one get for V 1 the profile equation 

(D + A)V 1 = 0, (268) 

which is a transport equation along the flow of v° with an extra term, local, of order 0. Notably 
this equation involves X only as a parameter. 

For the NS equations one get the profile equation: 

(D + A — ad 2 x )V l = (269) 

which is hyperbolic in t, x and parabolic in t, X. In particular the hypersurface {t = 0} is now 
non-characteristic . 

The equations (|268p - (|269p are therefore both well-posed when set for X in the whole real line. 
Let us explain an analogy which makes this sound natural. The expansions (|265j) - (|267|) could be 
seen as viscous and local counterparts of the asymptotic expansions of weakly nonlinear geometric 
optics (cf. for instance [93]). In this latter case the profiles are periodic with respect to the fast 
variable X and the small parameter \fv refers to short wavelengths. In this setting of geometric 
optics it is well-known -at least from experts- that a viscosity of size v only parabolizes the profile 
equations corresponding to an ansatz of the form (I265p - (|267p . without disturbing the well-posedness 
(on the contrary actually). We refer to the papers [63], [33], [M] which illustrate this remark with 
some close settings. 

Now if one looks for some solutions v v of the NS equations with a vortex patch vj as initial data, 
one has to prescribe zero initial data for the layers so that the compatibility condition between the 
transmission conditions and the initial condition on the "corner" {t = X = 0} are not satisfied even 
at order zero (cf. Remark l4.ip . This ruins any hope for smoothness with respect to X, which leads 
to some difficulties in the analysis of the stability of the expansions. One way to get the compatiblity 
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conditions is to choose the initial condition, restricting ourselves to a kind of well-prepared initial 
data, assuming that the initial data is already of the form (|265j) - ()267j) . that we could qualified by 
"well-prepared" initial data. It is possible to prove the existence of such well-prepared initial data 
thanks to some Borel lemma. We refer for an insight of the method to the papers |106j and [105] . 
For such data one observed that the transmission conditions persist when time proceeds, and the 
lifetime of such expansions is the one of the solution of the Euler equation ("the ground state") 
which traps the main part of the nonlinearity of the problem. 



4.11 Weaker singularities 

Let us mention that both kind of expansions (|200p - (|209p and f|265|) - ()267[) can be useful in the case 
of conormal singularities weaker than vortex patches. If for instance the vorticity is continuous 
through the initial internal boundary {(po(x) = 0} but not the derivative of the vorticity, then the 
profiles 0°, V 1 and P 2 do not depend on X anymore. The layers appear only at the following orders 
(as for instance in (|247jl ). More generally if for k G N the vorticity is C k through the initial internal 
boundary {ipo (x) = 0} it is possible to write a complete asymptotic expansion of the vorticity of 
the form: 

uj u (t,x)=uj (t,x)+Y^ t ) j n j (t,x) + yV^i 1 SP(t,x, y { ^ , ) + 0(VUF°). 

Such an observation was mentioned in the setting of the approximation of semi-linear symmetric 
hyperbolic systems of PDEs by the vanishing viscosity method in |105] . 



4.12 Stronger singularities 

Let us now talk about stronger singularities than vortex patches, for which layers of larger ampli- 
tude are expected. Let us first deal with vortex sheets which involve some discontinuity jumps of 
the velocity (instead of jumps of the vorticity). An initial velocity piecewise smooth with discon- 
tinuity jumps across an hypersurface {(fo = 0} is an extremely instable configuration for the Euler 
equations: in general the jumps of the velocity do not stay localized on a smooth hypersurface 
{<p°(t, .) = 0} when time proceeds. A few positive results are available with analytic initial data: 
local in time persistence was proved by Bardos and al. in |1L)8| (see also some extension to global 
persistence for small analytic perturbation by Caflisch and Orellana in [18J and by Duchon and 
Robert in [IH]), but the papers [T5], [HE], [H2]> [ID] ruin willingness to extend out the unphysical 
case of analytic data. Still a few physical phenomena are known to bring some stability such as 
compressibility (see [37] in the two dimensonal supersonic case) or surface tension (see [75], [8j). 

If for one of the previous reason the velocity v given by the Euler equation (or an appropri- 
ately modified inviscid system) stays piecewise smooth with discontinuity jumps only across an 
hypersurface {ip°(t, .) = 0}, we expect that the corresponding velocities given by the Navier-Stokes 
equation (or a modified viscous system) should be given by an expansion of the form: 

On ^ , tru _ <P°(t' x ) 



v u {t,x) ~ v"(t, x) + V(t,x, 



ut 



since the perturbation term V should have to compensate the jumps of the inviscid solution v° 
through the transmission conditions [V] = — [v ] where the first brackets denote the jump of V 
across {X = 0} and the second one the jump of v° across {(f°(t, .) =0}. In addition we have to 
prescribe [c?x^] = 0. Actually because of some nonlinear effects we have to consider also the next 
term in the expansion: 

v v (t,x) - v ^t,x) + V (t,x,^^) + V^iV 1 (t,x,^^) (270) 
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to close the profile equation. Actually plugging this ansatz into the NS equations and equalling 
according to the orders of y/ui yields -with the notations of the introduction- the equations: 

V° ■ n = 0, (L + -)V° = t{{V° ■ V x {v° + V°) + {V 1 ■ n)d x V } and dxV 1 ■ n = -div x V°. (271) 

This equation contains the same difficulty than the Prandtl equations: to get rid of V 1 in the second 
equation we use the third equation what leads to a loss of derivative. It is likely that here again 
some positive results are possible in the case of analytic data. Indeed Caflisch and Sammartino 
in [T7] have succeed to prove well-posedness for analytic data in the well-prepared counterpart of 
the equations (I27ip in two dimensions when the radius of curvature of the curve is much larger 
than the thickness of the layer. A solution more radical but also mathematically simpler to avoid 
this loss of derivative is to consider some cases without variation in the transverse directions. For 
instance we can consider the following very special case: we set the phase ip°(t,x) = x\ and an 
inviscid velocity v° of the form v°(t, x) = f(xi)e2 where the function / is on M-t. The solution 
of NS equations are then simply 

v »(t,x) = V(^=)e 2 , 
\Jvt 

where V is the solution of the elliptic equation 

X 
~2 



d\V + ^d x V = 



with the following transmission conditions on {X = 0}: 

[V] = -[/] and [dxV] = 0, (272) 

where [/] denotes the jump of / across {xi = 0}. 

If now we strengthen again the amplitude looking for expansions of the form 

\Jvt 

with a > 0, we get for V the equation 

X 
~2 



d\y + ^d x V + aV = 0, 



which is still coercive in H (M) for a < j but admit for eigenvalue for a = ^> with k G N*. The 

x 2 

corresponding eigenfunctions are the Hermite functions := d x ~ Hi with Hi := e~~ . Hence 
when strengthening the amplitude of the transition layer (that is when increasing a) there are 
still some non-trivial solutions of the profile equation but the nature of the profile problem totally 
changes, from a elliptic problem with non-homogeneous condition transmissions to an eigenvalue 
problem. 

In this latter case it is then possible to consider also some viscous perturbations singular with 
respect to several dimensions. For instance it is possible to exhibit some family of solutions of the 
two dimensonal NS equations with velocities of the form: 

v u (t,x) = -±=V(t,^=), (273) 

where the profile V(t, X) vanishes when the fast variable X := (Ai,X2) goes to infinity. The 
corresponding vorticities are of the form 

u u {t, x) = —Sl{t, -?L) with Q := curlx^. (274) 

Vl yjvt 
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Plugging these ansatz into the NS equations, whose vorticity formulation reads as follows in 2d: 

d t uf + v u ■ Vlj u = uAu v , (275) 
and equalling the terms of order jj^s and those of order we get the following pair of equations: 

V ■ V x fi = and (A* + ■ V x + 1)0 = tdtSl, (276) 

where we denote Ax := d x + d\ 2 and Vx := (dx 1 ,dx 2 )- We then see that the first equation in 
(|276p is satisfied if the vorticity profile £1 is radially symmetric since in this case the corresponding 
velocity profile V is orthoradial. The initial hypersurface {t = 0} is characteristic for the second 
equation in (1276j) which have therefore parasite solutions. Actually omitting the X dependence we 
get the ODE: 0, = tdfO, whose solutions are O = C\nt. However only one namely O = has a 
correct behaviour for t near 0. Now setting t = in the second equation we get the equation 

(Ax + \x • Vx + m = 

whose solution is given by Q(X) = Ce~ x2 ; where C is determined by the conservation of the 
vorticity mass. The ansatz (|273p describe the viscous smoothing of an initial Dirac mass at x = 
and are usually referred as Oseen vortex. We refer here for instance to the papers |54j . |53j . and to 
the references therein, for a much more precise study. We do not claim any novelty in the previous 
lines but we think that it was thought-making to incorporate a little bit of this material here. In 
particular we found interesting to gather the profile equation (12760 and the profile equation (|18j) 
corresponding to the vortex patches, and to observe in particular the shift of the spectrum caused 
by the difference of amplitude, which make the analysis quite different. 



5 Stability 

This section is devoted to the proof of Theorem 11.51 Once again we will work with the velocity 
formulation. We will prove: 

Theorem 5.1. For any k G N* there exists vq > such that for < v < vq for any (t,x) £ 
(0,T) x M. d 

v u (t,x)=v°(t,x) + yV^i J V j (t,x, ip °^ ) + V^i k+1 v u R 

U w 

where the profiles V 3 are the ones of the previous section and the sequence of the remainders 
(v R )o<u<u ^ in T. 

Theorem 11.51 is a corollary of this result. 
Proof. For any fc£N* there exists k' > k such that the function v v a given by 



:(t,x) := v°(t,x) + Y,V^t J V\t,x, fMl). 



3 

is an approximated solution of the NS equations in the sense that 



d t vl + v ■ VX + V xP » = uAv v a + Vrt k+1 F^ (277) 

div< = V^ fe+1 F£, (278) 
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with (FjQq <u<U0 and (F R )o <u<uo in T . We observe that to prove Theorem 15.11 it is sufficient to 

prove that v u R := y/ui k+1 (v u — v%) is in T. To do that we recast the NS equations into the 
following equations for the remainder v R : 

t{d t v R + v ■ V x v R ) + ^±±v R + V x p R = vtAv R + F R , (279) 

6ivv^ = F%, (280) 

where the NS pressure p v have been decomposed into 

\fui k+1 

p»(t, x) = p£(t, x) + *——jr R (t, x). 

We immediately express the pressure remainder p R in function of the velocity v R taking the diver- 
gence of the equation (|279p (and then using the equation (|280p ) so that the problem now reads 

t(d t v R + v ■ V x v R ) + ±±±v» R - vtAv R = n(/,^) + F R , (281) 

where II is bilinear of order and (F r )q <u<1/q is in T. It is now a classical observation since 
the paper [71] by Gues that -in such a process of remainder estimate- the nonlinear terms in the 
equation (I28ip have some vt in factor so that the (v R )o<i/<v live till T. As a consequence we now 
drop the tilde and the index v to simplify and we briefly explain how to get the needed energy 
estimate. We apply the operators Aj of spectral localization to the equations (128111 to get for 
j > — 1 the equations: 

t(d t + v ■ V x )A jVR + ^-A jVR - vtAAjVR = F RJ (282) 

where 

F RJ := AjU(v, v R ) + AjF R + t[v ■ V x , Aj]v R 

For p > 2 we multiply the equation (I282p by p| Ajt;^| p ~ 2 .Ajf^ and we integrate by parts in t, x on 
(0, T) x M. d to get the following estimate 



( {k + 1)p - 1) / |A^f + T( I \A jVR \ P ){T) +vtl P (p - lWAjvnflVxAjvnr 2 

1 J(0,T)xR d JR d J(0,T)xR d 



< / p\FR tj \.\A jVR \ P 1 <p||^j1li>((o,r)xR-)-l|Aj«flllL'( 1 fOT)xR'')(- 283 ) 

For s' G N, for any r' G (0,1), we multiply by 2-?( s ' +r >. Then we use the commutator estimates 
(I35|) and (|37H to control the first and the third terms of Frj. We take advantage of their factor 
t to absorb them by the second term in (|37p through Gronwall lemma. Finally we take the p-th 
squareroot and we let p goes to infinity to conclude that the sequence 

s'+r' 



v ~t ll^llL 00 ((0,T),C s '' r '(R 3 ))°< i/ < ;y o 

is bounded. Since this holds for arbitrary s' £ N and r' £ (0,1), this means that the family 
(v R )o<v<v is in T. □ 

Remark 5.1. Theorem \5.1\ proves the stability of complete expansions. It would be also interesting 
to look at the stability of short expansions that is to understand to what extent it is possible to 
justify an expansion with a few terms assuming as weak tangential smoothness as possible. In this 
direction another approach would be to develop a multi-phase extension of what is done here. 

I warmly thanks J-Y. Chemin, R. Danchin, T. Gallay and O. Glass for some useful discussions. 
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